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ABSTRACT 


Vibration  control  is  a  very  important  issue  in  satellites.  The  new  high-resolution  digital 
imaging  devices  are  especially  sensitive  to  vibrations.  Antennas  used  in  laser  communications  also 
require  a  very  quiet  environment  so  that  their  performance  is  not  degraded.  The  Stewart  platform 
is  capable  of  isolating  an  optical  payload  from  the  noisy  spacecraft  bus.  Until  recently,  only  pas¬ 
sive  methods  were  used  in  all  vibration  isolation  applications.  Recent  advances  in  Digital  Signal 
Processing  techniques  made  the  development  of  vibration  control  algorithms  possible,  but  these 
usually  require  large  computational  power.  This  work  explores  using  a  computationally  efficient 
vibration-isolation  method  for  optical  payloads  by  using  hexapods.  The  method  suppresses  the  vi¬ 
bration  at  the  assigned  frequencies  and  does  not  affect  unassigned  frequencies  if  the  plant  is  linear. 
The  mathematical  analysis  includes  convergence  analysis  and  the  effect  of  unassigned  frequencies 
in  the  output.  The  computational  requirements  of  the  algorithm  is  evaluated  and  is  compared  to  the 
Multiple-Error  Least  Mean  Square.  The  method  is  very  robust  to  nonlinearities;  its  performance  is 
comparable  to  the  Multiple-Error  Least  Mean  Square  with  a  fraction  of  the  computational  time  and 
memory  requirements.  Also  it  requires  very  little  plant  knowledge.  Theoretical  results  are  verified 
through  simulations  using  a  Single-Input/Single-Output  plant  and  a  nonlinear  hexapod  model.  The 
controller  was  also  experimentally  validated  in  two  different  hexapods  and  the  performance  was 
found  to  be  similar  to  or  better  than  the  performance  obtained  with  the  Multiple-Error  Least  Mean 
Square  method  when  a  noisy  reference  signal  is  used. 
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I.  INTRODUCTION 


A.  THE  VIBRATION  ISOLATION  PROBLEM 

Although  the  vibration  control  problem  has  been  studied  for  a  long  time,  it  received  increas¬ 
ing  attention  after  the  advent  of  submarines.  The  main  purpose  of  a  submarine  is  to  travel  unnoticed 
and  the  noise  generated  by  its  several  subsystems  can  defeat  its  purpose,  making  its  localization  and 
identification  easier.  As  a  result,  the  study  of  noise  and  vibration  received  much  attention,  especially 
during  the  Cold  War  when  several  techniques  for  noise  and  vibration  isolation  were  developed. 

During  the  sixties,  space  started  being  considered  as  a  very  strategic  area  and  imaging 
was  among  the  many  applications  that  found  their  way  to  the  military  spacecraft.  During  the  first 
decades  of  space  imaging,  vibration  was  not  as  important  as  it  is  today  because  of  the  limitations 
of  the  imaging  sensors  available  at  that  time.  Nowadays  satellites  are  larger  and  have  much  more 
rotating  machinery  than  in  recent  past.  Among  the  vibration-generating  devices  present  in  the  recent 
satellites  are  large  appendages,  which  snap  when  entering  or  emerging  from  an  eclipse,  solar  panels 
that  rotate  to  track  the  sun,  cryocoolers,  pumps,  reaction  wheels,  control  momentum  gyros  (CMG), 
reaction  jets  and  magnetic  torquers.  The  vibrations  produced  by  these  elements,  if  transmitted  to 
the  imaging  device,  can  degrade  its  performance  appreciably.  Figure  1 . 1  shows  an  example  of  such 
degradation. 

At  the  same  time,  the  electronic  industry  has  made  astonishing  advances  in  imaging  de¬ 
vices.  Resolution  is  now  much  higher  than  it  was  a  few  years  ago,  and  new  optical  techniques 
produce  lenses  and  mirrors  with  less  distortion,  generating  sharper  images.  Even  when  distortions 
are  present,  the  images  can  be  post-processed  on  the  ground,  using  the  latest  advances  in  image- 
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(a)  Jitter  free  (b)  Jitter  present  (Ifjrad) 

Figure  1.1.  Jitter  Effeet  on  an  Imaging  Deviee  (From  [1]). 


proeessing  teehniques  to  enhanee  them.  One  famous  example  is  the  spherieal  aberration  on  the 
Hubble  teleseope’s  original  primary  mirror.  Using  image  proeessing  teehniques,  the  teleseope  eould 
still  be  used  while  the  mirror  was  being  upgraded[6]. 

The  most  common  cause  for  vibration-induced  image  degradation  is  due  to  vibrations  pro¬ 
duced  elsewhere  and  transmitted  through  the  satellite  bus  to  the  imaging  device.  Isolating  the  noisy 
environment,  dampening  the  structure  or  isolating  the  instrument  itself  by  using  an  isolating  mount 
(see  Figure  1.2)  can  reduce  this  type  of  vibration.  The  vibration  can  also  be  produced  on  the  quiet 
side  when  moving  parts,  such  as  cryocoolers,  must  be  mounted  as  closely  as  possible  to  the  sensor. 

Passive  methods  are  used  in  most  vibration-control  and  vibration-isolation  implementations 
for  several  reasons:  they  are  robust;  they  do  not  consume  electrical  power;  they  are  very  reliable 
and  the  technology  is  mature.  Passive  approaches  can  be  applied  to  mounts  that  prevent  the  vibra- 
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Figure  1.2.  The  Vibration  Isolation  Problem  (From  [1]). 

tions  to  propagate  to  or  from  the  spaeeeraft  bus  or  to  inerease  the  dampening  of  a  larger  sttueture. 
Unfortunately,  passive  approaehes  do  not  solve  all  vibration  problems. 

When  a  vibration  isolation  mount  is  used,  the  vibration  has  two  main  paths  to  migrate  from 
the  noisy  to  the  quiet  side:  the  vibration  isolation  medium  and  the  umbilieal.  Even  though  it  is 
possible  to  design  very  soft  passive  mounts  to  preventing  the  vibration  to  propagate  through  the 
mount,  prevent  the  vibration  from  propagating  through  the  umbilieal  is  impossible.  In  the  same 
way,  passive  isolation  approaehes  are  ineapable  of  suppressing  the  vibration  generated  on  the  quiet 
side.  This  is  espeeially  important  for  infrared  imaging  for  whieh  the  sensor  must  be  eooled  and  a 
eryoeooler  must  be  plaeed  as  elosely  as  possible  to  the  sensor  in  order  to  inerease  its  performanee. 
Passive  solutions  ean  only  perform  vibration  eaneellation  in  eertain  eases. 

Disturbanee  isolation,  also  ealled  “vibration  suppression,”  is  the  suppression  of  periodie 
disturbanees.  These  disturbanees  are  usually  generated  by  rotating  maehinery  and  appear  as  lines 
in  the  frequeney  speetrum.  These  lines  have  a  power  eontent  mueh  higher  than  the  baekground 
vibration  noise.  If  the  frequeney  is  known  and  fixed,  TMDs  (tuned-mass  deviees)  ean  be  used,  but 
their  plaeement  on  the  strueture  is  eritieal.  Thus  TMDs  are  not  usually  suited  for  isolation  mounts. 
Aetive  methods  ean  deal  with  frequeneies  not  known  a  priory  or  with  frequeneies  fluetuating,  and 
the  plaeement  of  the  aetuators  is  not  as  eritieal  as  when  using  passive  methods. 
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B.  ACTIVE  VIBRATION  SUPPRESSION 


Since  meeting  the  requirements  using  passive-only  solutions  is  becoming  increasingly  dif¬ 
ficult,  active  solutions  have  been  pursued.  Current  research  indicates  that  the  most  viable  solution  is 
the  use  of  a  generic  local  vibration  isolator  that  combines  both  passive  and  active  isolation.  Among 
the  several  approaches  tested,  the  most  promising  is  the  Stewart  Platform  [7],  which  can  gener¬ 
ate  force  and  torque  in  any  direction.  It  also  exhibits  good  stiffness  properties  [2]  with  actuators 
much  smaller  than  those  needed  with  serial  manipulators.  The  need  for  localized  vibration  isolation 
prompted  the  development  of  the  Vibration  Isolation  and  Suppression  System  (VISS)  [8,  9,  3]  and 
the  Satellite  Ultraquiet  Isolation  Technology  Experiment  (SUITE)  [10]. 

A  hexapod  has  complex  dynamics  and  a  detailed  model  can  be  quite  complicated.  Even 
if  a  perfect  model  is  developed,  space  missions  usually  do  not  allow  the  opportunity  for  scheduled 
maintenance  in  order  to  repair  aging  or  failed  systems  and  thus  system  identification  may  be  needed. 
Eor  this  reason,  all  vibration  control  work  found  in  the  reviewed  literature  uses  some  sort  of  adaptive 
method.  This  approach  makes  it  possible  to  compensate  for  model  imperfections  and  aging. 

Since  the  hexapods  are  multiple-input/multiple-output  (MIMO)  systems  strongly  coupled 
[11],  expecting  that  the  best  performance  can  be  obtained  by  using  MIMO  controllers  is  reasonable. 
Some  methods  are  available  for  MIMO  plants,  such  as  Multiple-Error  EMS[12],  Clear  Box  [13,  4] 
and  Repetitive  Control  [14,  15].  Unfortunately,  the  computational  resources  available  in  satellites, 
limit  the  use  of  these  methods  for  vibration  isolation  of  optical  payloads. 

One  important  characteristic  common  to  all  hexapods  used  for  vibration  isolation  as  de¬ 
scribed  in  the  literature  is  the  use  of  smart  struts.  All  of  them  have  six  identical  struts  with  an 
actuator  and  a  vibration  sensor  embedded  (geophone  or  accelerometer),  meaning  that  using  one 
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independent  single-input/single-output  (SISO)  eontroller  for  eaeh  strut  is  possible.  This  approaeh, 


while  simple,  has  limitations  for  a  system  as  strongly  eoupled  as  a  hexapod,  and  mixed  results  have 
been  aehieved  [11,  16]. 

Eliminating  the  requirement  of  having  a  MIMO  eontroller  does  allow  more  ehoiees  sinee 
several  vibration-suppression  algorithms  were  never  extended  to  the  MIMO  ease.  In  2000,  Kuo  and 
Morgan  [17]  presented  an  interesting  review  of  some  of  the  vibration  isolation  algorithms  available 
and,  then  elassilied  them  into  “referenee-based”  and  “synthetie  referenee.”  The  first  group  uses  a 
referenee  signal  eorrelated  with  the  disturbanee  while  the  seeond  one  uses  some  other  information 
to  generate  an  internal  referenee  signal.  The  “synthetie  referenee”  group  eontains  is  one  set  of 
eontrollers  ealled  “noteh  filters.” 

Noteh  filters  have  several  properties  that  are  very  attraetive  for  vibration  isolation  in  satel¬ 
lites.  First,  no  external  signal  is  used  as  a  referenee.  Sinee  extra  information  about  the  frequeney 
is  given  or  generated  on-the-fly,  the  frequeney  information  ean  be  used  to  generate  a  synthetie  ref¬ 
erenee.  This  eliminates  the  requirement  of  having  an  additional  sensor,  at  the  expense  of  requiring 
additional  information  to  generate  the  referenee  signal.  Another  eharaeteristie  of  this  elass  of  eon¬ 
trollers  is  that  exeellent  results  ean  be  aehieved  for  suppressing  vibrations  generated  by  rotating 
maehinery,  also  ealled  tonal  disturbanees. 

Among  the  eontrollers  found  in  the  reviewed  literature,  the  one  used  in  1998  by  Bertran 
and  Montoro  in  [18]  (based  on  the  work  presented  by  Kuo  and  Min  in  1995  [19])  required  the  least 
amount  of  eomputational  effort  is 

This  eontroller  is  quite  simple  and  does  not  use  the  model  of  the  plant  explieitly.  Although 
this  eontroller  is  very  simple,  the  reviewed  literature  did  not  mention  this  eontroller  being  applied  to 
MIMO  systems,  espeeially  on  hexapods.  In  addition,  no  extension  of  Montoro  and  Bertran’s  work 
for  MIMO  plants  was  found.  More  important,  no  stability  analysis  was  found  in  the  literature  that 
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considers  a  generic  plant.  However,  such  analysis  is  essential  to  meet  the  high-reliahility  require¬ 


ments  of  space  applications.  Another  drawback  is  that  the  method  requires  the  knowledge  of  the 
frequency  to  be  suppressed.  Other  frequencies  could  be  present  on  the  error  signal  and  some  inter¬ 
ference  can  occur.  This  possibility  must  be  addressed  if  this  controller  is  to  be  implemented  in  an 
actual  space  system.  However,  this  problem  has  not  yet  been  addressed  in  the  literature  reviewed. 

Although  much  effort  has  been  spent  on  developing  vibration  control  and  isolation  methods, 
our  examination  of  the  literature  revealed  that  no  one  has  evaluated  the  computational  effort  that 
each  method  requires.  Consequently,  a  cost  analysis  must  be  developed  and  the  vibration-isolation 
method  selected  must  be  compared  to  the  Multiple-Error  LMS  implementation. 

Finally,  most  modeling  work  done  on  hexapods  was  from  a  robotics  perspective,  using 
direct  and  inverse  kinematics.  Li  and  Salcudean  presented  a  model  for  a  hexapod,  but  assumed  that 
the  dynamics  of  the  actuators  are  negligible  [20].  Lebret,  Liu  and  Lewis  derived  a  very  complete 
model  for  a  hexapod  starting  from  the  energy  approach,  but  unfortunately  the  geometry  assumed  in 
their  work  does  not  apply  to  the  hexapods  available  at  the  Satellite  Research  and  Design  Center: 
each  pair  of  actuators  was  connected  to  the  top  plate  through  a  single  joint.  In  1991,  Fujimoto, 
Kinoshita  et  al.  derived  a  model  from  the  robotics  perspective  [21],  but  assumed  that  “the  joint  does 
not  move  in  high  speed,”  preventing  its  application  from  being  used  for  vibration-isolation  purposes. 
A  complete  model  of  the  hexapod  suitable  for  simulations  is  needed  in  order  to  study  the  effects  of 
the  multiple-input/multiple-output  (MIMO)  nature  of  the  hexapod  without  the  interference  of  actual 
hardware  imperfections. 
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C.  PROBLEM  STATEMENT 


The  goal  of  this  research  is  to  study  a  computationally  efficient  algorithm  that  meets  the 
vibration-isolation  requirements  of  optical  payloads  using  hexapods.  The  method  needs  to  be  tol¬ 
erant  to  nonlinearities,  must  require  significantly  less  computational  resources  than  the  Multiple- 
Error  LMS  implementation,  and  the  method  must  perform  comparably  to  the  Multiple-Error  EMS 
method.  Additionally,  the  method  must  be  able  to  deal  with  multiple  frequencies. 

A  stability  analysis  and  a  study  of  the  effect  of  uncontrolled  frequencies  must  be  performed. 
The  results  must  be  validated  on  simulations  and  experimentally. 

Einally,  the  computational  cost  of  both  the  proposed  method  and  the  Multiple-Error  EMS 
approach  must  be  derived  and  must  be  compared  in  order  to  verify  if  the  proposed  method 
achieves  the  desired  goal  of  significantly  reducing  the  computational  cost  over  the  reference  method 
(Multiple-Error  EMS). 

D.  OVERVIEW 

Chapter  I  presents  the  vibration-isolation  problem,  reviews  the  literature  and  states  the  re¬ 
search  objectives. 

Chapter  II  describes  the  two  different  hexapods  used  to  perform  the  experiments. 

Chapter  III  provides  a  detailed  state-space  model  of  the  Precision  Pointing  Hexapod.  This 
model  allows  the  designer  to  derive  some  important  properties  of  the  plant  and  is  used  to  implement 
a  detailed  simulation  program. 

Chapter  IV  presents  the  Adaptive  Disturbance  Canceller  method,  derives  the  stability  cri¬ 
terion,  and  studies  the  method’s  performance  when  multiple  frequencies  are  present.  Simulations 
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are  used  to  validate  the  mathematieal  results  using  both  a  SISO  plant  and  the  model  developed  in 


Chapter  III. 

Chapter  V  derives  and  eompares  the  eomputational  requirements  for  both  the  Adaptive 
Disturbanee  Caneeller  and  the  Multiple-Error  LMS  approaeh,  as  well  as  derives  and  eompares  the 
eorresponding  memory  requirements. 

Chapter  VI  validates  the  results  experimentally  by  using  the  two  hexapods  available  at  the 
Satellite  Research  and  Design  Center,  at  the  Naval  Postgraduate  Sehool.  A  eomparison  between 
the  Adaptive  Disturbanee  Caneeller  and  Multiple-Error  EMS  is  also  presented. 

Einally,  Chapter  VII  presents  a  summary  of  the  work  and  eonelusions  and  some  suggestions 


to  overeome  the  main  weaknesses  of  the  Adaptive  Disturbanee  Caneeller. 


II.  HARDWARE  DESCRIPTION 


A.  OVERVIEW 

The  Stewart  platform,  also  known  as  a  “parallel  manipulator,”  was  introdueed  by  Stewart 
in  1965[7]  and  is  eomposed  of  two  plates  eonneeted  by  some  (passive  or  aetive)  links.  One  of  the 
plates  is  eonsidered  to  be  the  referenee  (base)  and  the  other  (top)  is  the  one  for  whieh  the  attitude 
and  position  are  to  be  eontrolled.  A  diagram  of  sueh  a  platform  is  shown  in  Figure  2.1.  The  platform 
ean  be  eontrolled  by  ehanging  the  length  of  the  links.  It  ean  be  shown  that  in  order  to  get  six  degrees 
of  freedom,  one  needs  at  least  six  aetuated  struts.  One  ean  also  verify  that  more  than  six  aetuators 
will  not  inerease  the  number  of  the  DOF  of  the  manipulator  (differently  from  the  serial-linked  ease). 
The  hexapods  have  several  interesting  eharaeteristies,  namely,  high  stiffness  and  redueed  size  and 
weight  when  eompared  to  serial-link  manipulators.  A  eomprehensive  study  of  hexapods  and  their 
eharaeteristies  ean  be  found  in  [2,  22,  23,  24,  25,  26]. 


Figure  2.1.  Hexapod  View  in  3D  (From[2]). 
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Section  B  describes  the  first  hexapod  available  at  the  Satellite  Research  and  Design  Center, 


the  Ultra-Quiet  Platform.  This  hexapod  is  equipped  with  high  bandwidth  piezoelectric  and  sensitive 
geophone  sensors  on  each  strut. 

Section  C  details  the  Precision  Pointing  Hexapod,  which  has  long  stroke  actuators  and  is 
capable  of  pointing  and  can  simultaneously  isolate  vibration.  Since  this  hexapod  is  modeled  in 
Chapter  III,  it  is  described  in  more  detail. 

B.  ULTRA-QUIET  PLATFORM 

This  Ultra-Quiet  Platform,  shown  in  Figure  2.2,  was  conceived  for  vibration-isolation  pur¬ 
poses.  It  follows  a  cubic  configuration  (see  Figure  2.3)  in  order  to  minimize  the  cross-coupling 
among  the  struts  [27].  The  hexapod  rests  on  a  satellite  bus  mockup  (Figure  2.4)  where  the  distur¬ 
bance  source  is  mounted.  The  whole  structure  is  mounted,  with  rubber  vibration  suppressors,  on  a 
Newport  RS4000  floating  table  in  order  to  provide  isolation  from  the  ambient  vibration  noise. 


Figure  2.2.  Ultra-Quiet  Platform. 
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Figure  2.3.  Cubic  Configuration  (From  [3]). 


The  Ultra-Quiet  Platform  struts  are  equipped  with  piezoceramic  stack  actuators  (PZT). 
These  actuators  have  a  very  large  bandwidth  (resonance  peak  around  \.2KHz)  and  a  displace¬ 
ment  of  ±50qm,  which  is  adequate  for  vibration  control.  The  power  to  the  actuators  is  provided 
by  a  PCB  Piezotronics  790A06  amplifier  (±200F  at  ±  100mA).  Each  smart  strut  is  also  equipped 
with  a  Geospace  GS-llD  geophone,  which  measures  linear  velocity.  These  sensors  have  a  natural 
frequency  of  l4Hz  with  a  damping  ratio  of  0.8.  The  signal  conditioning  is  provided  by  the  CSA 
Engineering  Active  Vibration  Control  System  (AVCS).  The  fixtures  have  passive  isolation  built-in, 
which  provides  damping  in  series  with  the  actuators. 

The  disturbance  source  produces  acceleration  parallel  to  the  z-axis,  as  defined  in  Figure  2.1. 
It  uses  an  Aura  Bass  Shaker  (AST-lB-4,  251T),  and  it  is  mounted  on  top  of  the  mockup,  below  the 
bottom  plate  of  the  hexapod.  The  shaker  signal  is  produced  by  a  DSP  board  and  is  amplified  by  a 
DC  power  amplifier  {Kepco  BOP  20-10M). 

The  implementation  of  the  algorithm  is  made  on  a  dSpace  Alpha  Combo.  This  system 
comprises  a  DS1103  and  an  Alpha  Board.  The  Alpha  Board  is  powered  by  a  DEC  Alpha  SQOMHz 
processor  with  2MB  of  RAM.  The  DS1003  board  acts  as  a  slave  I/O  board  and  contains  a  Texas 
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Instruments  TM320C40  processor  running  at  50MHz  and  5 12^^;  of  memory.  The  setup  is  shown  in 


Figure  2.4.  The  controller  runs  on  the  Alpha  Board,  while  all  the  I/O  functions  are  executed  on  the 
DS1003  board. 

All  the  I/O  are  performed  at  10^  samples / s  while  the  controller  and  disturbance  generator 
run  at  10^  samples /s.  Since  there  is  not  appreciable  power  contents  on  the  measurements  above 
2KHz,  anti-aliasing  filters  were  implemented  by  oversampling  the  measurements  (lOKHz),  passing 
them  through  a  Chebychev  (type  I,  corner  at  200Hz,  3  order)  and  finally  down-sampling  the  output 
of  the  filters  to  10^  samples /s. 

The  outputs  of  the  controller  are  up-sampled  to  10^  samples / s  and  then  low-pass  filtered 
by  a  type  I  Chebychev  filter  in  order  to  smooth  the  output  sent  to  the  actuators.  The  disturbance 
signal  received  the  same  treatment.  Filtering,  A/D  sampling,  D/A  conversions  and  communication 
between  the  controller  and  the  host  PC  occur  in  the  DS1003  board. 

The  code  is  implemented,  compiled  and  downloaded  to  both  boards  through  the 
dSpace/Matlab  Real-Time  Workshop  environment.  The  actual  controller,  as  well  as  the  disturbance 
generator,  is  custom-coded  in  C.  These  tasks  are  accomplished  by  using  a  host  PC,  an  Intel-based 
P-III  Dell  Dimension  XPS  600MHz-  The  same  host  PC  is  also  used  to  collect  all  the  relevant  data 
produced  during  the  experiment.  The  data  are  then  converted  into  Matlab’s  proprietary  format  for 
storage. 
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Figure  2.4.  UQP  Experimental  Setup  (From  [4]). 
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Figure  2.5.  Ultra-Quiet  Platform  Controller. 
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C.  PRECISION  POINTING  HEXAPOD 


The  main  limitation  of  the  Ultra-Quiet  Platform  is  its  small  stroke,  whieh  prevents  it  from 
being  used  for  steering.  In  order  to  evaluate  the  possibility  of  performing  both  pointing  and  vibration 
isolation,  a  hexapod  with  a  longer  stroke  is  needed.  To  address  this  need,  the  Satellite  Research  and 
Design  Center  aequired  a  new  hexapod:  the  Preeision  Pointing  Hexapod,  from  CSA  Engineering, 
Ine..  The  platform  is  designed  to  allow  both  vibration  eontrol  and  pointing.  In  order  to  aehieve 
the  large  stroke  required,  the  Preeision  Pointing  Hexapod  uses  voiee-eoil  aetuators,  whieh  provide 
±5mm  travel.  Aeeelerometers  are  mounted  inline  to  provide  vibration  information.  The  platform, 
as  delivered  by  CSA  Engineering,  Ine.,  is  depieted  in  Eigure  2.6. 

This  hexapod  was  seleeted  for  modeling  beeause  it  has  a  larger  stroke.  This  hexapod  will  be 
used  for  pointing  researeh  as  for  well  as  vibration-isolation  researeh  and  thus  an  aeeurate  model  is 
needed.  A  eomplete  state-state  model  is  an  invaluable  tool  for  testing  new  algorithms  before  imple¬ 
menting  them  on  the  aetual  hardware.  In  order  to  have  the  model  developed,  a  deep  understanding 
of  the  hardware  is  needed  and  therefore  the  deseription  of  the  Preeision  Pointing  Hexapod  hardware 
is  more  detailed  than  the  understanding  needed  to  simply  operate  it. 

1.  Actuators 

This  platform  is  well  suited  for  position  and  vibration  eontrol.  The  aetuators,  manufaetured 
by  Moltran  Industries,  Ine.  (model  AFX  70N),  have  a  stroke  of  more  than  ±5mm  and  ean  provide 
more  than  2.5^  of  tilting  and  10"  of  twisting.  The  aetuators  ean  deliver  up  to  40A  of  dynamie 
foree  and  up  to  70A  of  statie  foree,  being  adequate  to  eontrol  large  levels  of  vibration.  The  most 
important  eharaeteristies,  from  Moltran’s  datasheet,  are  shown  on  Table  2.1  and  a  pieture  of  one  of 
the  aetuators,  as  mounted,  ean  be  seen  in  Eigure  2.7. 
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Value 

Unite 

Model 

Axial  Force  Transducer  AFX  VON 

- 

Force  af  dOWoc,  coil  af  25'’C 

100 

N 

Force  af  40Wrms,sme,  coil  af  25"C 

70 

N 

Force  af  40Wdc?  coil  af  140^C 

70 

N 

Force  af  40Wrms,sme,  coil  af  140^C 

50 

N 

Gain  af  fhe  Origin 

17.5 

N 

A 

Gain  af  ±5mm 

15 

N 

A 

Table  2.1.  Actuator  Specifications. 


The  actuators  are  powered  by  a  custom  power  supply  provided  by  CSA  Engineering,  Inc, 
composed  by  six  independent  switching  inverting  amplifiers  (swifching  frequency  of  29KHz).  The 
amplifier  has  negafive  gain  (—1)  and  accepfs  a  volfage  in  fhe  range  ±4E  on  fhe  inpufs  fhrough 
coaxial  cables.  The  oufpufs  of  fhe  individual  amplifiers  are  direcfly  connecfed  fo  fhe  acfuafors. 

2.  Accelerometers 

The  Precision  Poinfing  Hexapod  is  equipped  wifh  six  acceleromefers  {Kistler,  8304B2 
K-Beam)  mounfed  inline  wifh  fhe  sfruf’s  axis,  as  seen  in  Figure  2.1.  These  acceleromefers  have 
a  range  of  ±2g  and  are  linear  over  fhe  range  0  fo  200 //z.  The  main  specifications  are  shown  in 
Table  2.2. 


3.  Disturbance  Generator 


A  vibrafion  source  is  provided  in  order  fo  simulafe  periodic  disfurbances.  This  is  accom¬ 
plished  by  using  an  Aura  Bass  Shaker  AST-2B-4  mounfed  on  a  custom  adapfer.  The  shaker  charac- 
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Value  Unit 

Range  ±2  g 

Sensitivity  1000  mV  /g 

Frequeney  Response  ±5%  300  Hz 

Resolution  0. 1  mgrms 

Table  2.2.  Accelerometer  Specifications. 

teristics  are  shown  in  Table  2.3.  The  shaker  adapter,  shown  in  Figure  2.8,  allows  vibrations  along 
the  z-axis,  x-y  axis  (as  defined  in  Figure  2. 1),  fwisf,  fip  and  combinafions  of  fhe  previous.  The  shaker 
can  also  be  mounfed  on  fhe  lop  plale,  producing  vibralion  along  fhe  z-axis. 

The  signal  lhal  drives  fhe  shaker  is  produced  by  fhe  DSP  conlroller  and  senl  fo  an  AC  power 
amplifier  Sony  STR-D  Amplifier.  The  oulpul  of  fhe  amplifier  is  connecled  fo  fhe  shaker. 

4.  Floating  Table  and  Environment  Isolation 

The  whole  platform  is  sofl-mounled  on  a  Newporl  Slabilizer^^  laminar  flow  isolalor,  pre- 
venfing  fhe  ground  vibralion  from  interfering  wilh  fhe  experimenl.  In  order  fo  maximize  fhe  shaker’s 


Figure  2.8.  Base  Shaker  Adapter. 
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Value 

Unit 

Model 

AST-2B-4 

Magnet  Type 

eeramie 

Power  Rating 

50 

W 

Foree,  Nominal,  at  Resonanee 

132 

N 

Weight 

1.125 

Kg 

Resonanee  Frequeney 

42 

Hz 

Frequeney  Range 

20-100 

Hz 

Table  2.3.  Shaker  Speeifieations. 


effieieney,  the  lower  plate  is  mounted  on  the  NewPort  table  using  very  soft  rubber  pads.  Sinee  these 
pads  flatten  due  to  the  hexapod’s  weight,  they  are  replaeed  before  eaeh  set  of  experiments  to  ensure 
the  highest  degree  of  isolation. 

5.  Electronic  Support 

Both  the  power  supply  and  the  anti-aliasing  filters  are  implemented  in  the  “Control  Box.” 
This  box  aeeumulates  the  funetions  of  power  supply  for  sensors,  shook  deteetor,  emergeney  shut¬ 
down,  signal  eonditioning  and  anti-aliasing  filter  and  noise  suppressor.  Figure  2.9  shows  the  main 
funetions  implemented  in  the  Control  Box. 

Power  is  provided  by  a  eustom  off-the-shelf  power  supply  that  provides  ±5F  and  ±12F, 
whieh  are  the  voltages  that  both  the  box’s  eleetronies  and  the  sensors  need. 

The  aeeelerometers  are  eonneeted  by  a  standard  DB25  eonneetor  that  provides  them  with 
-I-12F  and  ground.  It  also  reeeives  the  aeeelerometer’s  signals  («  W /g),  whieh  are  sent  to  the 
anti-aliasing  filters. 
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fc/50 


Figure  2.9.  Control  Box  Main  Functions. 

A  sample  from  the  signal  from  the  accelerometer  #1  is  sent  to  a  window  comparator.  The 
normal  operation  of  the  hexapod  should  not  exceed  ±2g,  which  is  the  range  of  the  accelerometers.  If 
the  top  plate  of  the  hexapod  exceeds  ±2g,  the  accelerometer’s  output  voltage  will  be  either  smaller 
than  0.5F  or  larger  than  4.5F  and  thus  a  shock  is  considered  to  have  occurred.  In  response  to  this 
condition,  a  dedicated  circuit  will  disconnect  the  actuator  and  disturbance  amplifiers  from  the  1  lOV 
rail,  cutting  their  power  supply.  The  power  can  only  be  restored  by  pressing  the  reset  button.  The 
same  effect  can  be  obtained  by  pressing  the  panic  button.  This  mechanism  is  essential  to  preserve 
the  hexapod  integrity  if  the  controller  becomes  unstable,  especially  at  lower  frequencies. 

The  anti-aliasing  filters  are  implemented  as  fourth-order  switched-capacitor  Butterworth 
filters.  The  corner  frequency  can  be  adjusted  by  changing  the  frequency  of  the  oscillator  by  means 
of  a  potentiometer.  The  corner  frequency  is  defined  as  fc  =  fosc/50.  Since  fhis  is  acfually  a  dig- 
ifal  fiber  (alfhough  if  is  implemenfed  wifh  analog  componenfs),  one  musf  ensure  fhaf  fhere  are  no 
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frequencies  above  /c/2  when  selecting  fosc-  The  cutoff  frequency  was  selected  at  50Hz  due  to  the 
strong  nonlinearities.  The  filtered  measurement  signals  are  then  sent  to  a  connection  panel  where 
lOuF  capacitors  are  used  to  reduce  the  noise  even  further.  Those  capacitors  are  needed  to  prevent 
the  environment  noise  from  corrupting  the  measurements. 

6.  Controller  and  Host  Computer 

The  hexapod  is  controlled  by  a  dedicated  PowerPC-based  board  {dSpace  DS1103),  which 
is  installed  inside  the  host  PC.  A  connection  panel  is  attached  to  the  board  and  provides  all  the 
external  communication  with  the  board.  The  relevant  interactions  are  shown  in  Figure  2.10. 

The  host  computer  is  a  PC  Dell  Dimension  XPS  T500.  The  controller  is  implemented  by 
using  dSpace  and  Mathworks’  Real  Time  Workshop.  All  filtering,  disturbance  generation  and  I/O 
operations  are  implemented  with  Simulink  block  diagrams.  The  actual  control  laws  are  implemented 
in  either  C  or  Simulink.  All  the  data  logging  is  performed  using  dSpace"^  Control  Desk,  which  al¬ 
lows  the  desired  variables  to  be  streamed  to  disk  and  then  converted  to  Matlab’s  proprietary  format. 
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Figure  2.10.  Precision  Pointing  Hexapod  Experiment  Setup. 


III.  MATHEMATICAL  MODEL 


A.  OVERVIEW 

This  chapter  presents  the  development  of  a  mathematical  model  for  the  Precision  Pointing 
Hexapod  suitable  for  vibration  simulation  and  control.  This  model  enables  testing  of  different  con¬ 
trol  laws  without  having  to  implement  them  in  the  actual  hardware,  and  thus  enables  the  engineer 
to  evaluate  the  actual  influence  of  each  parameter  on  the  controller  performance. 

•  Section  C  lists  the  symbols  used  in  the  derivation. 

•  Section  D  derives  the  equations  used  by  the  model  starting  from  the  momentum  conservation 
equations,  according  to  the  Newtonian  formulation. 

•  Section  E  describes  the  Euler  angles  and  Euler  Parameters,  explaining  why  the  Euler  angles 
were  selected  for  the  hexapod’s  model. 

•  Section  E  modifies  the  free-body  equations  to  take  into  account  the  fact  that  the  forces  acting 
on  the  top  plate  are  assumed  to  act  only  on  the  joints. 

•  Section  G  incorporates  the  actual  actuator’s  model. 

•  Section  H  provides  the  measurement  equations  and  thus  completes  the  model. 

B.  INTRODUCTION 

Controlling  a  system  requires  knowledge  of  its  behavior.  A  control  system  may  even  operate 
without  a  priori  knowledge  of  the  system,  but  it  must  build  this  information  over  time  to  perform 
successfully,  as  in  the  adaptive  case.  Even  for  adaptive  controllers,  a  good  model  of  the  plant  is 
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important  in  order  to  provide  adequate  parametrization.  Another  reason  for  a  mathematieal  model 


is  given  by  the  nature  of  the  measurements,  sinee  they  are  usually  not  the  quantities  direetly  involved 
in  the  requirements.  In  faet,  the  desired  measurements  for  the  platform  are  roll,  piteh,  yaw  and  the 
angular  and  linear  aeeelerations.  None  of  these  quantities  are  measured  direetly  by  the  sensors:  they 
are  related  to  the  sensors  by  a  mapping,  whieh  ean  be  statie  or  dynamie. 

For  the  partieular  problem  of  the  hexapod  subjeet  of  this  researeh,  linear  and  angular  ae¬ 
eelerations  ean  be  statieally  derived  from  the  aeeeleration  of  eaeh  joint  of  the  moving  plate  as  a 
linear  transformation.  However,  angular  positions  (|),  9  and  \|/  (roll,  piteh  and  yaw)  are  related  to  the 
aeeelerations  by  differential  equations.  To  provide  estimates  of  the  desired  variables,  a  model-based 
teehnique,  sueh  as  the  Kalman  Filter,  has  to  be  implemented.  The  performanee  and  reliability  of  its 
results  are  as  good  as  the  reliability  of  the  mathematieal  model  of  the  system. 

The  rest  of  this  seetion  deseribes  in  detail  all  the  steps  involved  in  the  model  eonstruetion 
and  presents  the  results.  A  detailed  deseription  of  how  the  equations  were  derived  is  very  important 
so  that  the  same  methodology  ean  be  applied  to  derive  models  of  other  hexapods. 

C.  LIST  OF  SYMBOLS,  NOMENCLATURE  AND  DIMENSIONS 

The  Preeision  Pointing  Hexapod  was  presented  in  Seetion  II.C  and  only  relevant  details 
are  ineluded  in  this  seetion.  The  basie  geometry  of  the  platform  is  shown  in  Figures  3.1  and  3.2.  A 
hexapod  is  geometrieally  defined  by  its  12  joints  (six  upper  joints  and  six  lower  joints).  The  origin  is 
defined  as  fhe  geomefrie  eenfer  of  fhe  six  upper  joinfs  when  fhe  plafform  is  af  resf  wifh  fhe  aefuafors 
af  mid-eourse  (Figure  3.1).  The  axes  of  fhe  eoordinafed  system  are  also  depleted  in  Figure  3.1. 

The  physieal  quanfifies  used  fo  define  fhe  sysfem  in  fhe  formulalion  are  shown  in  Table  3.1. 
The  definilion  of  fhe  symbols  used  in  fhe  derivafion  are  in  Table  3.2. 
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Figure  3.2.  Hexapod — Top  View. 
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Symbol  Value 


.25 


r  0.0277 

0.0 

0.0 

0.0 

0.0277 

0.0 

L  0.0 

0.0 

0.0553 

100271 

Kg 

Kgm^ 


0.1271 

104271 


Description 


Radius  of  the  circle  that  contains  all 
the  upper  joints  (see  Fig.  3.2) 

Radius  of  the  circle  that  contains  all 
the  lower  joints  (see  Fig.  3.2) 

Angular  separation  of  the  clos¬ 
est  lower  joints,  measured  from 
the  center  of  the  lower  joints  (see 
Fig.  3.2) 

Angular  separation  of  the  clos¬ 
est  upper  joints,  measured  from 
the  center  of  the  upper  joints  (see 
Fig.  3.2) 

Distance  between  upper  and  lower 
planes,  when  the  platform  is  at 
“zero”  (see  Fig.  3.1) 

Mass  of  the  top  plate 

The  moment  of  inertia  of  the  plat¬ 
form 

Natural  frequency  of  the  actuator 


[none]  Damping  factor  of  the  actuator 

[none]  Gain  factor  for  the  input  to  force  re¬ 
lation  in  the  actuator 

rad/s  First  pole  of  the  accelerometer 

rad/s  Second  (and  last)  pole  of  the  ac¬ 
celerometer 
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Table  3.2.  Symbols  Used  in  the  Mathematieal  Model. 
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D.  FREE-BODY  MODEL 


In  order  to  make  the  problem  more  traetable,  it  is  assumed  that  all  forees  and  torques  are 
measured  with  respeet  to  an  inertial  referenees  system  (defined  in  Seetion  C).  The  most  basie 
equations  for  a  free -body  model  are  the  equations  of  momentum  eonservation  (Equations  (3.1) 
and  (3.2)).  These  qualities  do  not  ehange  if  no  external  torque  or  foree  is  applied  to  the  body.  The 
linear  (Q)  and  angular  (H)  momentum  are  defined  as: 


Q  =  m\ 


(3.1) 


H  =  Jw, 


(3.2) 


where  m  is  fhe  mass,  J  fhe  inertia  matrix,  w  the  angular  veloeity  and  v  the  linear  veloeity  (all  with 
respeet  to  the  top  plate). 

The  foree  F  and  torque  T  aeting  on  the  top  plate  are 

F  =  Q  =  m\  +  m\  (3.3) 

T  =  H  =  jw  +  Jw.  (3.4) 


Equations  (3.3)  and  (3.4)  ean  be  simplified  if  the  mass  and  the  inertia  matrix  ean  be  eon- 
sidered  eonstant,  as  it  is  the  ease  when  the  equipment  mounted  on  the  top  plate  ean  be  eonsidered  a 
rigid  body.  With  this  assumption  rh  =  0  and  j  =  0  and  Equations  (3.3)  and  (3.4)  ean  be  simplified 
to 


my 

(3.5) 

Jw. 

(3.6) 
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E.  EULER  ANGLES 


There  are  two  main  approaehes  that  ean  be  used  to  define  the  angular  veloeity:  Euler  angles 
and  Euler  parameters  (also  known  as  quaternions). 

The  first  approaeh  is  adequate  for  angles  within  bounds,  sinee  this  approaeh  has  singularity 
problems.  Eurthermore,  it  uses  transeendental  funetions  (sin,  eos  and  tan),  making  its  eomputation 
expensive.  Using  the  azimuth-elevation-roll  sequenee  and  aeeording  to  [28],  the  derivatives  of  the 
Euler  angles  are 


1 

tan  9  sin  (|) 

tan  9  eos  (|) 

0 

= 

0 

eos  9 

—  sin  9 

¥ 

0 

sin(l) 

COS(|) 

COS0 

COS0 

w  =  M(a)  w 


(3.7) 


and  the  singularity  is  at  9  =  7t/2.  The  platform  movements  are  small,  with  all  angles  being  less  than 
15°,  and  the  singularity  will  never  be  reaehed.  Without  the  singularity,  the  only  drawbaek  of  the 
Euler  angles  approaeh  is  its  eomputational  requirements. 

The  other  possibility  is  to  use  Euler  parameters  (quaternions).  Aeeording  to  [28],  they  ean 
be  written  as 


e'l 

Wxei  +  Wj.ea  +  Wze4 

e'l 

1 

-Wxei  -  -  Wye4 

e's 

“  ~2 

-Wyei  -f  w^e2  -  Wxe4 

e'A- 

-Wj-ei  -  Wye2  -  Wxe-i 

with  the  eonstraining  equation 


^1 +^2  +  ^3  +  ^4  =  1’ 


(3.9) 


where  e\,  e2>  ^3  and  ca,  are  the  quaternions. 
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This  approach  has  two  main  advantages  over  using  Euler  angles:  there  is  no  singularity  and 
the  updates  require  only  standard  multiplications.  The  Euler  angles  approach  was  selected  for  this 
model. 

Assuming  that  the  body  and  inertial  frames  have  the  same  origin,  a  vector  in  the  inertial 
frame  °p  can  be  transformed  into  the  body  frame  using  the  transformation  [28] 

c9c(|)  —5^9 

=  —  celts'll/ +  5(|)59c\|/  c(|)c\|/ +  i'(|)i'95\|/  5(|)c9  °P  (3.10) 

+  c(|)59c\|/  —  i'(|)c\|/  +  c(|)595\|/  c(|)c9 

=  '’R/P,  (3.11) 


were  c  stands  for  cos  and  5  for  sin.  The  transformation  matrix  is  orthonormal  and  invertible  for 
the  small  angles  that  the  Precision  Pointing  Hexapod  can  produce.  The  inverse  transformation  is 

Using  the  Euler  angles  and  the  assumption  of  constant  mass  and  inertia  matrix,  a  simple 


model  can  be  written  as 


V 

m 

tan  9  sin  (|)  tan  9  cos  (|) 


cos  9 

—  sin  9 

sin(l) 

COS(|) 

COS0 

COS0 

J“^T 


(3.12) 
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where  the  foree  F  and  the  torque  T  aet  on  the  eenter  of  mass  of  the  top  plate. 

This  model  deseribes  a  body  floating  in  free  spaee  and  has  arbitrary  foree  and  torque  as 
inputs.  Therefore,  this  system  is  eontrollable  as  is  its  linearized  version.  The  main  problem  with 
this  model  is,  however,  that  the  foree  and  torque  are  not  eontrolled  direetly:  they  are  generated  by 
the  six  aetuators. 


F.  FORCES  ACTING  ON  THE  TOP  PLATE 


Disregarding  the  gravitational  effeets,  there  are  only  six  external  forees  and  no  external 
torque  aeting  on  the  top  plate  (through  the  joints).  A  mathematieal  relationship  must  be  found  that 
ean  translate  these  forees  into  the  resultant  foree  and  torque.  It  is  assumed  that  the  aetuators  are 
mass-less  and  the  forees  generated  by  the  aetuators  are  parallel  to  the  aetuator’s  axis  and  aet  only  at 
the  joints. 

The  foree  aeting  on  the  eenter  of  mass  is  the  summation  of  all  forees  aeting  on  the  body. 
Defining  L,  as  the  foree  aeting  on  the  joint  i  and  e  =  [ei  as  the  foree  produeed  by  the 

aetuators,  the  resultant  foree  is 


F  =  XC 


i=\ 


f  b,-B, 


B, 


(3.13) 


where  e,  is  the  foree  produeed  by  the  aetuator  i,  b,  is  the  position  of  upper  joint  i  and  B,  is  the 
position  of  the  lower  joint  i. 


b,  =  p  +  r, 


=  p  +  ^Ror^,,- 


(3.14) 
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Defining  iT,  =  the  foree  F  ean  be  written  as 

F  =  id 

i=i 

I 


7x  7i  53  54  55  56 


=  5(p,a)e. 


(3.15) 


The  torque  ean  be  evaluated  in  a  similar  manner.  Eaeh  foree  has  a  assoeiated  torque  when 
translated  to  the  eenter  of  mass  and  these  torques  ean  then  be  summed  to  produee  the  total  torque 
T.  First,  define  a  veetor  r,-  as 


r, 


Rr,r„ 


(3.16) 


where  p  is  the  position  of  the  eenter  of  mass  of  the  top  plate  and  b ,  is  the  position  of  the  joint  i. 


Therefore,  the  overall  torque  ean  be  written  as 


T 


b,  -B,  \ 
b,  -B,  \ 


'F(p,a)  e. 


(3.17) 


where  5,  =  r,-  x 


b,-Bi 

[bi^ 
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The  results  of  Equations  (3.15)  and  (3.17)  provide  a  linear  relationship  between  the  input 


forees  and  the  resulting  foree  and  torque.  By  using  these  results,  the  plant  model  ean  be  rewritten  as 

P 

V 

d 

w 

where  the  parameters  of  iT  (p,OC)  and  OT  (p,a)  were  dropped  for  eompaetness  and  are  implieit. 

G.  ACTUATORS 

The  foree  of  the  aetuators  eannot  be  eontrolled  direetly.  The  eontroller’s  output  is  a  num¬ 
ber  that  is  eonverted  to  a  voltage/eurrent  and  then  sent  to  the  aetuator.  The  aetual  foree  produeed 
depends  on  the  aetuator’s  internal  dynamies. 

The  aetuators  installed  on  the  Preeision  Pointing  Hexapod  are  voiee-eoils,  whieh  are  ba- 
sieally  a  eoil  in  the  inside  and  a  magnet  on  the  outside  (Figure  3.3).  A  eurrent  flowing  in  the  eoil 
generates  a  magnetie  field  that  interaets  with  the  magnetie  element  ereating  the  foree,  whieh  is  linear 
for  small  displaeements.  Both  a  spring  foree  and  a  viseous  foree  must  also  be  taken  into  aeeount. 
Therefore,  the  equation  relating  the  foree  generated  by  the  aetuator  to  the  eurrent  flowing  through  it 
and  the  position  is 

F{t)  =  Kii{t)  —  Kxx{t)  —  bx{l)  —  max{t) ,  (3.19) 

where  Ki  is  the  gain  from  eurrent  to  foree  and  b,  nia  and  Kx  are  the  viseous  frietion,  moving  mass 
and  spring  eonstant  of  the  aetuator. 


M(a)  w 

J-^Te 


(3.18) 
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Ax  I  s 


Figure  3.3.  Actuator  Diagram. 


The  force  produced  by  the  actuator  must  be  written  in  terms  of  the  inertial  coordinate  sys¬ 
tem.  The  displacement  x  is  the  distance  traveled  by  the  joint  along  the  unity  vector  ||{^|~g'||  =  jT; 

Xi  =  J'i 

=  (p  +  r,-bo,;)^5; 

=  (^p  +  ^Roro,,  -bo,;)  J,-. 

Since  b^,;  =  ro,;,  this  becomes 

^  =  (p  +  "Robo,-bo,,y5,-.  (3.20) 

The  velocity  v,  of  the  joint  i  is  the  sum  of  the  velocity  of  the  center  of  mass  and  the  velocity 
of  the  joint  with  respect  to  the  center  of  mass.  Since  the  latter  is  produced  by  the  angular  velocity 
of  the  top  plate,  the  resulting  expression  is 

Vi  =  (v  +  wxr^j^jTi 

=  (\  +  yfx‘’R„ro,t^  Ji.  (3.21) 
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Similarly,  the  acceleration  of  the  joint  i  is  the  acceleration  of  the  center  of  mass  summed  to 


the  centripetal  acceleration: 

a,-  =  a  +  w  X  (w  X  r,-) .  (3.22) 

Projecting  a,  it  on  the  struts’  axes  leads  to: 

Ui  =  (a  +  wx  (wxr,))^jr,- 

Ui  =  (^a  +  wx  ^wx^R^ro,;^^  iT,-.  (3.23) 

Substituting  x,  x  and  x  into  Equation  (3.19),  the  force  on  joint  i  is 

Fi{t)  =  Kiii{t)  -  +  ^Robo,,'  -  bo,;^  +b(\  +  y/x  + 

(^a  +  w  X  X  ^  ^  5,-,  (3.24) 

assuming  it  is  parallel  to  the  actuator’s  axis.  All  the  terms  of  this  equation  are  known  states  except 
for  the  input  /,(t). 

It  is  important  to  mention  that  the  output  of  the  controller  is  not  a  current,  but  a  voltage.  This 
voltage  is  sent  to  a  DC  amplifier  and  its  output  goes  to  the  actuators.  The  amplifier’s  bandwidfh  is 
200Hz  and  fhe  fime  consfanf  of  fhe  acfuafor’s  LR  circuif  is  Im^.  Therefore,  fhe  model  developed  in 
Ibis  secfion  is  valid  only  up  fo  200Hz-  If  higher  frequencies  are  needed,  fhen  fhe  elecfrical  dynamics 
musf  also  be  modeled. 

1.  Model  With  Actuators 

The  force  Equafion  (3.24)  can  be  used  in  Equafions  (3.12)  fo  generate  a  more  complete 
model.  To  do  fhaf,  Equafion  (3.24)  musf  be  wriffen  in  a  more  compacf  form.  Defining 

hi  =  (^Kjc  +  ^Robo,,'  -  bo,;^  +b(\  +  v/x  + 

+  w  X  X  ^RoCo,,-^  )  )  ^'’ 


35 


h^[hl 


(3.26) 


and 


u=  [/i  ...  kf . 


Equation  (3.12)  can  be  rewritten  as 


p 

V 

V 

m'^ 

a 

M(a)  w 

w 

r^'Te 

Fi...Fe 

M(a)  w 


Fi...Fe 


iyWi'-h) 

M{a)w 
[j-i‘r(^,u-h) 
with  h  and  u  as  defined  in  Equations  (3.25),  (3.26)  and  (3.27). 


(3.27) 


(3.28) 


H.  ACCELEROMETERS 


All  hexapods  used  for  vibration  control  reviewed  in  the  literature  have  some  kind  of  sensor 
to  measure  the  vibration  built  into  the  strut.  The  Precision  Pointing  Hexapod  is  not  an  exception 
and  has  accelerometers  built  inline  with  the  actuator’s  axis.  Section  ll.C.l  (p.  15)  describes  the 
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mounting  bracket  and  the  accelerometers,  while  Table  3.3  shows  the  relevant  specifications  of  the 
accelerometers. 


Value 

Unif 

Range 

±2 

g 

Sensifivify 

1,000 

mV/g 

Resonanf  Freq.  (nom.) 

1.2 

kHz 

Freq.  Response  (±5%) 

300 

Hz 

Finearify 

±1 

%FS 

Transverse  Sensifivify  (nom.) 

<  1 

% 

Phase  Shiff  (lOHz,  nom.) 

2 

degree 

Resolution  (0.5  fo  lOOHz) 

0.1 

rugrms 

Table  3.3.  Accelerometer  Specifications. 


As  seen  in  Section  G,  the  model  is  only  valid  up  to  200Hz  due  to  restrictions  of  the  actuator’s 
model.  Since  200Hz  is  within  the  accelerometer’s  bandwidth,  their  dynamic  can  be  considered  just 
a  constant.  It  is  also  important  to  realize  that  the  accelerometers  mounted  on  the  struts  measure  DC 
level,  and  this  model  does  not  limit  this  behavior. 

The  accelerometers  measure  the  acceleration  of  the  joint  projected  on  the  actuator’s  axis. 
The  measurement  can  then  be  written  as 


at  =  ^  9 i 

b;  -B; 


(3.29) 


where  a,  is  the  measurement  of  accelerometer  i  and  a,  is  the  acceleration  at  the  joint  i.  This  acceler¬ 
ation  was  evaluated  in  Equation  (3.23). 

Using  the  above  result  and  defining  as  fhe  acceleromefers’  consfanf  gain,  fhe  measure- 
menf  equations  become 


m 


ai 

at. 

Ka,  (a  +  wx  (wx*Roro,i))^iri 


Ka,  (a  +  w  X  (w  X  *Roro,6) )  ^  ^6 


(3.30) 
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Therefore,  the  last  equation  is  the  measurement  equation  for  the  state-spaee  model  de- 
seribed  in  Equation  (3.28),  whieh  ean  be  used  for  vibration  simulation  purposes. 
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IV.  VIBRATION  SUPPRESSION 


A.  CHAPTER  LAYOUT 

This  chapter  contains  the  details  of  the  Adaptive  Disturbance  Canceller.  First,  Section  B 
derives  the  Multiple-Error  LMS  method. 

Section  C  presents  an  overview  of  the  Adaptive  Disturbance  Canceller.  The  mathematical 
details,  including  a  convergence  analysis,  are  presented  in  Section  D. 

Note  that  the  Adaptive  Disturbance  Canceller  isolates  the  vibrations  at  the  assigned  fre¬ 
quencies  only.  As  a  result.  Section  D.3  presents  a  study  of  the  behavior  of  the  controller  in  the 
presence  of  unassigned  frequencies  on  the  error  signal. 

Section  E  presents  simulations  results  obtained  with  the  Adaptive  Disturbance  Canceller. 
The  model  derived  in  Chapter  III  was  used  to  simulate  the  PPH  in  Section  E.2.  The  Adaptive 
Disturbance  Canceller  is  a  SISO  controller  and  multiple  controllers  might  interact  with  one  another 
on  a  distributed  controller  configuration.  Thus,  full  nonlinear  model  simulations  are  of  special 
interest. 

B.  MULTIPLE-ERROR  LEAST  MEAN  SQUARED 

Advances  in  Digital  Signal  Processing  techniques  led  to  several  estimation  techniques  and 
many  of  these  are  also  valuable  for  noise  cancellation.  Techniques  commonly  applied  in  noise 
cancellation  are  usually  based  on  adaptive  algorithms.  The  general  approach  is  shown  in  Eigure  4.1. 

Among  adaptive  filter  implementations  used  to  cancel  the  disturbance  signal  d  [n]  one  of  the 
simplest  and  best  understood  is  the  Eeast  Mean  Square  method  (see  [29]  for  a  detailed  description 
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Figure  4.1.  Noise  Caneellation  Problem. 


and  speeifie  properties).  This  filter  assumes  an  all-pole  plant  {H)  and  minimizes  the  eost  funetion 

7  =  E{|eHn,  (4.1) 

where  E{«}  is  the  expeetation  operator.  Additionally,  the  LMS  is  a  first  order  method,  whieh 
eontributes  to  its  simplieity  and  redueed  eomputational  effort  in  eomparison  with  more  sophistieated 
algorithms. 


Figure  4.2.  An  Adaptive  Vibration  Controller. 


The  vibration  problem  has  one  additional  eomplieation  when  eompared  to  the  noise  ean- 
eellation  problem:  the  transfer  funetion  H  is  after  the  filter.  In  faet,  the  filter  provides  the  inputs  to 
the  plant  H.  An  updated  diagram,  refieeting  the  vibration  eontrol  problem  is  shown  in  Figure  4.2. 
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It  is  important  to  realize  that  Figure  4.1  and  Figure  4.2  only  differ  on  the  position  of  the  plant  and 


the  eontroller:  they  are  switehed.  Unfortunately,  the  very  faet  of  switehing  the  plant  and  the  filter 
prevents  a  straight  LMS  from  being  used  for  eontrol.  The  LMS  filter  has  a  weight  update  equation 
in  the  form 


W[n+l]  =  W[n\+iue*[n\u[n\,  (4.2) 

and  sinee  the  plant  is  plaeed  after  the  filter,  a  ehange  on  u  [n]  will  not  refleet  on  e  [n],  but  on  e  [u  +  1] 
(the  LMS  assumes  an  all-pole  plant).  Therefore,  e  [n]  does  not  depend  on  u  [n]  and  the  LMS  method 
eannot  be  applied  direetly  to  the  vibration  eontrol  problem. 

Widrow[30]  and  Burgess[31]  solved  this  problem  working  independently  in  1981.  Their 
method  used  an  estimated  model  of  the  plant  and  is  ealled  Filtered-X  LMS.  This  algorithm  inherits 
most  of  the  eharaeteristies  of  the  simple  LMS  scheme.  Unfortunately,  the  Filtered-X  LMS  method 
does  not  address  one  important  issue:  MIMO  systems.  Elliot,  Stothers  and  Nelson  generalized  the 
original  Filtered-X  LMS  in  1987  [12],  presenting  the  Multiple-Error  EMS.  The  new  algorithm  was 
able  to  control  MIMO  plants,  but  still  possesses  the  restriction  of  having  a  single  signal  correlated 
with  all  disturbances.  In  essence,  the  Multiple-Error  EMS  transforms  the  MIMO  plant  into  a  SIMO 
plant,  with  a  single  input  {x{n])  and  several  outputs  {eo{n])  (see  derivation  in  Section  1).  The 
structure  of  the  controller  proposed  in  [12]  can  be  seen  in  Eigure  4.3. 

Both  the  Multiple-Error  EMS  and  the  Eiltered-X  EMS  approaches,  as  well  as  their  variants, 
seem  to  be  the  most  common  active  vibration  controFisolation  algorithms  in  use  today.  They  are 
among  the  least  expensive  in  terms  of  computational  requirements  and  are  comparatively  simple  to 
implement.  The  Multiple-Error  EMS,  due  to  its  widespread  acceptance,  is  used  in  this  work  as  a 
benchmark  and  is  presented  next. 
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Figure  4.3.  LMS  for  Vibration  Control. 


Figure  4.4.  Multiple-Error  LMS. 


1.  Derivation  of  the  Multiple-Error  LMS 


The  Multiple-Error  LMS  assumes  that  the  MIMO  plant  H  (Eigure  4.4)  ean  be  modeled  as 
a  Einite  Impulse  Response  (EIR)  filter  or  order  M  with  Q  outputs  and  P  inputs.  Thus,  the  q*^  output 
[n\  due  to  the  input  p  is  given  by: 


M 

yq,p  W]  =  L  H  “P  [n-m], 

m=0 


(4.3) 
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assuming  that  the  Q  disturbances  dq  [n]  acting  on  the  system  are  correlated  to  the  signal  ;c  \n] ,  which 
must  be  measurable. 

Further,  assuming  W  to  be  a  matrix  of  dimensions  S  x  P,  the  signal  u  [n]  can  be  generated 
as 

5-1 

“p  [n]  =  [”]  (4.4) 

.v=0 

Substituting  the  input  value  derived  in  Equation  (4.4)  into  Equation  (4.3),  the  output  of  the 
plant  at  time  n  is 

P-l  M 

yq  VA  =  L  L  ^q^p  C'”]  “p  “  '"I 

p=0m=0 

P-l  M  5-1 

=  12  12  Vp  H  12  ^■'^P [n-m]x[n-m-s\ 

p=0m=0  A— 0 

P-l  M  5-1 

=  12  12  Y^Kp[p\ws,p['^-m\x[n-m-s].  (4.5) 

p=0m=0.s— 0 

The  error  can  be  then  computed  as 


Cq  [n]  =  dq  [n]  -  Jq  [n] 

p-l  M  5-1 

=  <5?^  [«]  -  XI  12  12  ^q^P  w  ^*:P  [n-m\x[n-m-s\.  (4.6) 

p=0m=0 .9=0 

Eollowing  the  same  reasoning  that  is  used  in  the  derivation  of  the  EMS  method,  W  must  be 
updated  minimizing  the  cost  function 


7  =  E 


e-i 

12  4  ['A 


q=0 


(4.7) 
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By  using  the  steepest  deseent  and  as  the  learning  rate,  the  weight  ean  be  adjusted  as 


follows: 


Ws,pW\  =  Ws,p[n]-ii 


^Ws,p  [n] 


=  Ws,p[n]-^i- 


'dws,p  [n] 


Q^d^{e\[n]} 

=  ws,p[n\ -  A'  L  ^ — nr' 

„=0  ^Ws,pVl] 


Unfortunately,  eomputing  the  expeetation  E  [n\ }  in  real  time  is  impossible,  but  the  same 
approximation  used  in  the  derivation  of  the  standard  LMS  algorithm  ean  be  made  here: 


[n]. 


By  using  the  above  approximation  and  proeeeding  with  the  ealeulation  of  Ws^p  [n  +  1] : 


w,,p  [n  +  1]  «  Ws,p  [«]  -  A/  2,  ^  r 

q=0  ^^s,p[n\ 

^  r  1  o  r  1 

Pi  w,^p[n\-2^l^eq[n\^^ 

q=0  OWs,p  [« 


(4.10) 


Next,  substituting  the  error  Cq  [n]  from  Equation  (4.6)  into  Equation  (4.10),  Wg^p  [n-\- 1]  ean 


be  evaluated  as 


Ws,p  [«  +  1]  ~  Ws,p  [n]  + 


.  v.e-1  r 


(4.11) 
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Assuming  that  the  weights  ehange  slowly  (p  small),  then,  Ws^p  [n  —  m]  ^  Ws^p  [n]  and  the 


previous  equation  ean  be  simplified  to 


Ws,p  [n+l]  «  Ws,p  [n]  + 


hLj  /[«— mkf/t— m— 

^2-1  r  1  V^/y=0^«J=0 v=0  q,p''-  ^  J  L  ij 

Ws,p  [n]  + 


.  ^0-1  r  1  V  /^=0  5^=0  L 

2^11.%  e,[n]  ^ - 


Ws,p  [n]  + 


o  y'G"!  r  1  ( v/’— 1  vM  f  r  1 

2A'I:J.0  ^qH  (^ly=0  ^-=0  ^y=0  Kp'  [H 


x\  n—m—s 


(4.12) 


Note  that  the  derivative 


tion  (4.12)  eollapses  to 


Sw'./.p'H  ; 


'dw,[p[n\ 


is  different  from  zero  only  if  p  =  p'  and  s  =  s'.  Therefore,  Equa- 


e-i  M 


Ws,p  [u  +  1]  ~  Ws,p  [n]  -f-lpYj^q  W]  Yj  ^‘i’P  [m]x[n-m-s]. 


(4.13) 


Now,  defining  rq^p  [n]  as 


'q,p 


fhen,  Equafion  (4.13)  beeomes 


^=0  m=0 


M 

[n]  =  Y,  hq^p[m]x[n-m] 
m=0 


e-i 


(4.14) 


Ws,p  [n+\]Ki  Ws,p  [n]  +2pY^q  ^q^p  [n-s], 

q=0 


(4.15) 


Equafion  (4.15)  is  fhen  used  fo  define  fhe  weighf  updafe  equafion  for  fhe  Multiple-Error 


EMS: 


e-i 

Ws,p  [n+f]  =  Ws,p  [n]  +2pY‘^q  M  ''q,p  [n-s]. 

q=Q 


(4.16) 


One  imporfanf  eharaeferisfie  of  fhe  Multiple-Error  EMS  mefhod  is  fhaf  a  planf  model  is 
needed  fo  generafe  Pq^p  [n].  Aefually,  fhe  ferm  Y,m=o^q,p  [p]x[n  —  m  —  5]  ean  be  inferprefed  as  being 
fhe  q'^  oufpuf  of  fhe  planf  when  fhe  p'^'  inpuf  is  x  [u  —  5]  and  all  fhe  ofher  inpufs  are  zero. 
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C.  PROPOSED  METHOD 


For  a  time-invariant  or  slowly  varying  tonal  disturbance,  having  a  full  model  of  the  plant 
or  a  controller  designed  to  work  over  a  frequency  band  is  an  overkill.  The  system  is  only  excited 
at  a  particular  frequency  and  thus,  the  transfer  functions  are  reduced  to  single  complex  numbers. 
Assuming  that  the  plant  is  linear  and  stable,  it  is  possible  to  generate  any  arbitrary  sinusoidal  at  the 
output  of  the  system  and  thus  cancel  any  sinusoidal  disturbance.  A  controller  using  this  approach 
was  presented  in  1998  by  Bertran  and  Montoro  [18].  The  proposed  controller  needs,  for  the  studied 
case,  only  a  few  assumptions:  a  stable  linear  SISO  plant,  tonal  disturbances  with  known  frequencies 
and  the  plant  without  zeros  at  the  frequencies  of  interest.  A  block  diagram  describing  the  controller 
is  shown  in  Figure  4.5. 

The  controller  proposed  can  be  summarized  as  follows:  assuming  that  the  plant  H  is  linear 
(see  Figure  4.5),  then,  for  any  sinusoidal  signal  d[n]  with  frequency  (Oc  it  is  possible  to  find  a 
sinusoidal  input  x  [n]  such  that  y  [n]  =  —dn  (if  (e^®^)  ^  0).  Usually  one  would  write  the  input  as 

x[n]  =  AcOs((Oc?^  +  Pa:)-  (4.17) 

Several  algorithms  which  compute  the  optimal  value  of  x„  by  minimizing  e[n\=y  [n]  +  d  [n] 
are  available,  but  almost  all  assume  that  y  [n]  is  a  linear  combination  of  the  parameters.  The  input 
Equation  (4.17)  can  be  changed  to  the  equivalent  form 

x[n]  =  acos((Oc«) +  ^sin(c0c«)-  (4.18) 

Next,  assuming  that  d{  (e^®"^)  =  ae^P,  the  steady-state  output  y  [n]  can  be  written  as 

y[n]  =  aacos(c0c?^  + P) +^ctsin(c0c?^  + P)-  (4-19) 

Using  this  form,  the  output  y[n]  is  linear  in  the  parameters  a  and  b  and  most  adaptive 
algorithms  can  be  used  to  find  a  and  b.  Berfran  and  Monforo  used  an  updafe  similar  fo  fhe  one 
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Figure  4.5.  Adaptive  Disturbanee  Caneeller. 


defined  in  the  LMS  method: 


a[n-\-\]  =  n[n] +/re[n]eos(cOc?^)  (4.20) 

b[n-\-\\  =  b[n\ -\- jJLe[n\?,m{G)cn) ,  (4.21) 

with 

e[n]=y  [n]  +  d[n].  (4.22) 

D.  MATHEMATICAL  ANALYSIS 

1.  Stability 

Bertran  and  Montoro  derived  the  stability  analysis  for  the  setup  presented  in  this  paper,  but 
the  partieular  plant  assumed  is  not  appropriate  for  the  Stewart  Platform  problem.  Therefore,  a  more 
general  approaeh  is  presented  in  the  following  diseussion. 

The  disturbanee  d  [n]  ean  be  written  as 

d[n]  =AtlF(e-^'“^)e-^'“^”,  (4.23) 
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where  cOc  is  the  frequeney  of  the  disturbanee. 

Assuming  q  small,  A  [u]  ehanges  slowly  and  the  output  y  [n]  ean  be  written  as 

y  [u]  A  [n]  ^  (e^'“^)  (4.24) 

Thus,  sinee  e[n\=y  [n\  -\-d[n],  the  error  is 

e\n]  =  d  \n]  +  y  \n\ 

=  (A  +  A[u])l//(e-^'®-)e^'“^” 

=  A[u]l/f(e-^'®^)e-^'®^".  (4.25) 

Using  the  weight  update  equation  and  the  definition  of  A  [u]  =  A + A  [u] ,  the  update  equation 
for  A  [n\  is  given  by 

A  [u  +  1]  =  A  [u]  +  qe  [n]  e“-^®‘'”. 

Adding  A  to  both  sides  of  the  previous  equation  leads  to 

A  [u  +  1]  =  A  [u]  +  qe  \n]  e“-^®‘'”.  (4.26) 

Also  note  that  algorithm  eonvergenee  may  be  proven  if  |A  [u  +  1]  |  <  |A  [u]  |,  for  all  integer 
values  of  n.  Therefore, 

|A[u+1]|^  =  A*  [u+ l]A[u+ 1] 

=  (A*  [n]  +  qe*  [n]  e-^'®^")  (A  [n]  +  qe  [n]  e--^'®^") 

=  I A  [n]  I V  A*  [u]qe  [n]  e”-/'®^"  +  A  [n]  e^'®^”  +  \^\^\e[nf .  (4.27) 

Now,  using  the  error  value  from  Equation  4.25,  the  above  expression  is  redueed  to 

I A  [u  +  1]  I  ^  =  |A  [u]  I  ^  +  A*  [n]  liA  [n]  d{  (e^®^)  e^'®-"e“^'®^”  + 

A  [nYA*  [n]  di*  (e^®^)  e-^'®^”e-^'®^"  + 

Iql^A*  [n]  d{*  (e-^®^)  e--^'®^"A  [n]d{{t^^^)  e-^'®^". 
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Applying  complex  conjugate  properties,  the  final  expression  for  |A  [n+  1]  ^  is  given  by 


I A  [n  +  1]  I '  =  \A[n]f  +  lRe{A  [nYA*  [n]  (e^'“^) }  +  \^\^  \A[n]f  \!H  (e^'“^ )  | ' 

=  |A[n]|^(l  +  2/?e{A/:f/(e-/'“^)}  +  H^|^(e-''“0|^)’  (4-28) 

and  it  can  be  concluded  that  lim„_><x>  |  A  [n]  |  ^  =  0  provided  that 

1  +  IRe  (e^'“^) }  +  \n\^  (e^'“^)  <  1, 

or,  equivalently, 

2Re  {fi9{  (e-'®^) }  +  |qp  \9{ (e-^®^)  <  0. 

The  weight  A  [n]  converges  to  —A,  as  A  [n]  =  A  +  A  [n] .  Thus,  the  error  converges  to  zero. 
Assuming  that  (e^®^)  |  <<  1,  the  stability  criterion  can  be  simplified  to 

Re  {qlff  (e^®^) }  <  0.  (4.29) 

As  a  guideline,  q  can  be  selected  by  using 

9{*  (e-/®-) 

(4.30) 

10|lf/(e^“0r 

where  the  magnitude  of  the  plant  is  not  critical  and  the  phase  must  be  known  within  ±7t/2.  This 
guideline  is  not  strict  and  the  model  of  the  plant  does  not  need  to  be  precisely  known. 


2.  Optimal  Weights 


Assuming  that  the  controller  converges,  then  y  [n]  =  —d  [n].  Assuming  that  the  disturbance 

d[n]  is 


r/[n]  =D,e-^'(®^"+^4, 


(4.31) 
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then,  the  output  y  \n]  can  be  evaluated  as 


e[n]  =  y  [n]  +  d  [n] 


The  convergence  of  the  algorithm  ensures  that  the  error  converges  to  zero  and  A  \n]  con¬ 
verges  to  a  constant.  Therefore, 


e[n\  =  y  [n]  +  d  [n] 


0  =  Al//(e-^®^)e^®^"  +  D, 


j&A  J&ctr  I  p.  j{iacn+yc) 


Defining  d{  (e-^®^)  =  ae^P,  the  previous  expression  becomes 


=  +  aAe^'(®^”+P) 


n  e^'(®^”+T^-)  =  -aA*e^'(®^"+P) 


(4.32) 


where  A*  is  the  optimal  weight. 

The  real  coefficients  a*  and  are  obtained  by  rewriting  Equation  (4.32)  in  the  real  form 

a*  =  -^cos(Yc-P) 


b*  = 


■sin(Yc-p). 


(4.33) 


3.  Crosstalk 


It  is  important  to  verify  how  the  algorithm  behaves  in  the  presence  of  uncontrolled  fre¬ 
quencies.  For  such  an  analysis  the  disturbance  is  assumed  to  contain  two  complex  exponentials  at 
frequencies  00^  and  (0„c-  Thus 


d  [n]  =  Dce-^'®^"+^'^^  + 


(4.34) 
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where  cOc  is  the  disturbanee  to  be  suppressed  and  co„c  is  not  to  be  suppressed. 


It  is  assumed  that  €  N  (where  N  denotes  the  set  of  Natural  numbers)  sueh  that  A  [n]  = 
A[n  +  N]  and 


<  ytg  4  e  and 


The  input  of  the  plant  ean  be  written  in  the  eomplex  exponential  form  as 


x[n]  =A[n]e^'“^". 


(4.35) 


Taking  the  DFT  of  x  [n]  and  applying  the  frequeney  shift  property,  X  {k)  beeomes 


X{k)  =  AL{k-kc). 


The  output  y  \n]  is 


y[n]  =  X \n]  *h\n]. 


and  therefore,  {y  [«]}  is 


‘y{k)  = 

=  Al{k-kc)9{(&^^^^y 


(4.36) 


The  weights  A  [n\  are  updated  as 


A[n+ 1]  =  A[n] +qe[n]e 


(4.37) 
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Since  e{n\=y  [n]  -\-d{n],  Equation  (4.37)  can  be  rewritten  as 

A  [n  +  1]  =  A[n\-\-  jjL{y[n]-\-  d[n]) 

=  A[n]+  lue-j^^’^y [n]  +  +  D,e-/'®-"+^'^-) 

=  A[n]  +  lue-j^^’^y  [n]  +  qD,e^'^^e-''®^'’e-^'“^”  +  qD„ce-/'^-e^'®-"e--''®^” 

=  A  [n] +  qe“'^®‘'”3^[n] +qDce-^"''‘' +qD„ce-^"''"^e-^^“”.  (4.38) 

Applying  the  DFT  to  Equation  4.38,  leads  to 

®)F‘r{A(n+ 1)}  =  ®)r‘r{A[n]}+q!D)r‘r  {e“-^'®^”y[n]}  + 

©ITT  { 1 }  +  qD„ce^'^"^  ©ITT  |  e-^'^“”  | . 

Using  DET  properties,  the  previous  equation  becomes 

•  iKk 

j^{k)eJ^  =  AL{k)+iu‘X{k  +  kc)+ijDceJ^^5{k)+ijDnce^'^'^5{k-ks).  (4.39) 

By  substituting  “X  (k)  from  Equation  (4.36)  into  Equation  (4.39)  leads  to 
AL  (k)  =  J4  {k)  +  {k)  d{  ^  5  (^)  +  5{k  —  k^).  (4.40) 

The  final  expression  can  then  be  obtained  by  grouping  the  terms  in  AL  (k) 

AL{k)  (ej^  -I- ^e^'(T+®c)^^  =  qDce^'^^5(A:)  + 

ljD„,&^^'-5{k-k8).  (4.41) 

Equation  (4.41)  shows  that  the  weights  have  components  only  at  ^  =  0  and  k  =  k^  and 
thus  the  error  only  has  components  at  (0^  and  (0„c-  Therefore,  the  controller  does  not  create  new 
frequencies.  Influence  of  the  controller  in  the  error  at  co  =  co„c  is  evaluated  by  replacing  ^  =  ^5  in 
Equation  (4.41)  which  leads  to 

Al  {k^)  ^e^“^  —  1  —  /jdf  ^e^(  "  ^  ^  =  qDce^^'"5 (^g)  +  tiDnc&^'^''‘'^{k^  —  kf) 

=  (4.42) 
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Thus, 


mh)  = 


jjnt 


Qj-fr  —  1  —  iji!}{ 


Similar  results  may  be  obtained  by  evaluating  Equation  (4.41)  at  =  0 


(4.43) 


^  (0)  (e^'^  -  1  -  qi/C  ^e^'(  ^+“^) )  )  =  qD^e^'^'^S  (0)  +  qD„ce-^'^"^5  (0  -  h) 

=  (4.44) 


Thus, 


MO) 


(4.45) 


Using  the  faet  that  ‘X{k)  =  J^{k-  kc)  ^  •E{k)  =  7' {k)  +  ‘D {k)  ean  be  evaluated 


as 


£  (fc)  =  7  {k)+  DnM'^'^  5  (^  -  ^s)  +  DcM'^  5  {k) 

=  J^{k  —  kc)  ^  +  Dnc&^'^'^5{k  —  knc)  +  DcQ^^"?> (k  —  kc) 


(  D 

I  — J—. — -b{k  —  kc)  + 


DncQ^'^''‘^^{k  -  knc)  +  DcM‘^5  {k  -  kc) . 


•  27C*s 

qJ^t  —  1  —  qj//'(e-/“"^) 


+ 


(4.46) 


At  this  point.  Equation  (4.46)  shows  that  the  error  has  eomponents  at  COc  and  (One  only. 
Eurthermore,  the  disturbanee  at  (One  does  not  affect  the  error  at  (0^-  Accordingly,  the  disturbance 
at  cOc  does  not  affect  the  error  at  the  frequency  (One-  Thus,  the  controller  does  not  generate  new 
frequencies. 

Evaluating  Equation  (4.46)  atk  =  knc  the  error  due  to  the  noncontrolled  disturbance  is 

‘Eiknc)  =  .2.,,  - i;/(et®-)+D„,e^'^-.  (4.47) 

qJ  n  _  1  _^j;/'(ei®nc) 
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Assuming  that  So  is  small  {k^  «  N),  then  e  "  1  and  Equation  (4.47)  can  be  simplified 


to 


£  {knc)  = 


Jjnc 


iKk^ 


qJ  N  —  1  —  (ei®nf  ) 


(e-^®”)  +  D„ce^'''"‘ 


— qff/ (e-t®"f) 

f  1  - 


1 


9{{ej^nc) 


^  (e-^“"") 


0. 


(4.48) 


Equation  (4.48)  shows  that  the  error  tends  to  zero  when  the  assigned  frequency  tends  to 
the  exact  disturbance  frequency  (6co  tends  to  zero).  This  important  result  shows  that  there  is  not 
discontinuity  in  the  solution  in  the  vicinity  of  the  true  solution,  provided  (0„c  is  not  a  zero  or  pole  of 
the  plant. 

However,  assuming  6q)  is  sufficiently  large  and  (e^®"^)  |  <<  1,  leads  to 


nT) 

^(k  \  -  ^  _ 

^yi^nc)  —  ,2Ttk^ 

qJ^t  —  1  —  (ef®«f) 


^  +  Dnce^^''‘ 


fjDfic^ 


Jlnc 


e-'  A"  —  1 


(e-^®"^)  + 


1  +  ^  I  . 


(4.49) 


e-^  «  —  1 

So,  using  (e-^®"^)  |  <<  1  and  assuming  that  ||e^^“  —  1||  1,  Equation  (4.49)  may  be 

approximated  by 


£  {knc)  ~  1  + 


pi  f^jjnc 

i^nr^ 


.2Kk^ 


H  (e-^®"") 


(4.50) 


The  above  result  shows  that  the  controller  has  very  little  impact  on  the  component  at  (0„ 
when  this  frequency  is  far  from  00^  and  q  is  small. 
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E.  SIMULATIONS 


Simulations  are  conducted  to  verify  the  mathematical  development  and  to  verify  whether 
the  assumptions  made  are  reasonable.  Recall  that  stability  conditions  were  derived  assuming  a 
single  frequency  on  the  error  signal,  and  thus  it  is  important  to  verify  the  behavior  of  the  controller 
if  other  frequencies  are  present  in  e  [n].  Additionally,  the  Adaptive  Disturbance  Canceller  is  a  SISO 
controller;  therefore,  it  is  necessary  to  verify  if  it  can  control  a  hexapod,  as  it  is  a  strongly  coupled 
MIMO  plant. 

1.  Single-Input/Single-Output  System 

The  first  simulation  is  intended  to  evaluate  the  algorithm  convergence  when  the  distur¬ 
bance  consisted  of  a  single  tone.  The  method  was  implemented  in  Simulink  as  described  by  Equa¬ 
tions  (4.20)  and  (4.21).  The  plant  used  was  a  second  order  system  with  a  natural  frequency  of  100//z 
and  damping  of  3%  and  its  output  was  sampled  1,000  times  per  second.  The  controller  also  ran  at  a 
\KHz  sampling  rate.  White  noise  was  added  to  the  error  signal,  with  SNR  =  20dB.  The  results  are 
shown  in  Fig.  4.6. 

Results  show  that  the  controller  works  as  expected,  even  for  such  a  high  noise  level.  In 
addition,  note  that  the  controller  does  not  affect  the  system  at  uncontrolled  frequencies. 

The  next  simulation  studies  the  controller  behavior  when  two  tones  are  present  in  the  dis¬ 
turbance  signal  (fc  =  33Hz  and  /„c  =  AlHz),  where  33Hz  is  the  disturbance  to  be  suppressed.  Fig¬ 
ure  4.7  shows  that  the  controller  cancels  the  assigned  frequency  and  does  not  change  the  unassigned 
one  {AlHz)-  The  weights  exhibits  components  at  /  =  QHz  (from  fc),  f  =  \AHz  and  /  =  %QHz-  The 
latter  two  values  are  predicted  by  theoretical  results,  since  \fd  —  fc\  =  l^Hz  and  fc  +  fd  =  80//z. 
Figure  4.8  shows  the  performance  of  the  algorithm  when  noise  is  added  to  the  sensors  {e  [n]). 

Theoretical  results  predicted  that  this  controller  cannot  be  applied  directly  when  the  phase 
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f(Hz) 

—  controller  on  - controller  off 

Figure  4.6.  Single  Noisy  Tone. 

of  the  plant  is  |  +  Ikii  because  1 1  +  /rae-'P  |  =  1 1  +  jpa  sin(  p)  |  >  1 .  In  order  to  verify  how  sensitive  to 
this  parameter  the  method  is,  the  plant  was  changed  so  that  its  phase  at  fc  is  almost  90°.  Figures  4.9 
and  4.10  show  the  resulting  performance  of  the  controller.  As  expected,  the  maximum  |^|  is  very 
small  and  thus  the  convergence  rate  is  very  slow. 

The  last  SISO  simulation  tests  the  performance  of  the  algorithm  when  several  controlled  and 
some  uncontrolled  frequencies  are  present.  Figure  4. 1 1  demonstrates  the  results  for  four  controlled 
(10,  14,  34,  and  36Hz)  and  three  uncontrolled  (13,  35  and  40Hz)  tones.  Noise  was  added  to  the  error 
signal.  It  is  relevant  to  mention  that  among  the  controlled  frequencies  there  were  34,  35  and  36Hz, 
where  only  34  and  36Hz  are  to  be  suppressed.  The  goal  was  to  verify  how  selective  the  algorithm 
is.  Results  show  the  tone  at  35Hz  is  left  untouched  while  the  tones  at  34  and  36Hz  are  attenuated. 
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Figure  4.7.  Case  B,  No  Noise. 
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Figure  4.8.  Case  B,  Noisy. 
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Figure  4.9.  Case  C,  No  Noise. 
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Figure  4.10.  Case  C,  Noisy. 
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—  controller  on  - controller  off 

Figure  4.11.  Case  D:  Four  Controlled  and  Three  Uncontrolled  Tones. 
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2.  Multiple-Input/Multiple-Output  System:  The  PPH 


The  Precision  Pointing  Hexapod  is  a  highly  coupled  system  as  it  is  impossible  to  move  any 
node  separately  from  the  others.  The  Adaptive  Disturbance  Canceller  was  developed  as  a  SISO 
controller  and  either  the  controller,  must  be  generalized  to  a  MIMO  approach,  or  a  multiple-SISO 
approach  must  be  taken  (also  called  decentralized  control).  Previous  attempts  in  using  decentralized 
control  had  mixed  success.  Beavers  [11]  did  not  obtain  a  good  performance  for  higher  frequencies 
and  the  performance  for  lower  frequencies  was  poorer  than  the  one  obtained  by  using  the  MIMO 
methods  (Multiple-Error  LMS,  Repetitive  Control  and  Clear  Box  [4,  14,  13]).  Good  results  using 
distributed  control  (six  independent  SISO  controllers)  were  obtained  with  some  controllers  [27,  32]. 
Since  previous  work  exhibited  mixed  results,  the  mathematical  model  developed  in  Chapter  III 
(p.  23)  is  followed  to  verify  if  the  coupling  is  a  limiting  factor  for  the  use  of  the  Adaptive  Disturbance 
Canceller  on  hexapods. 

The  simulation  setup  was  exactly  the  same  as  that  used  in  Section  E  (p.  55),  but  a  multiple- 
SISO  implementation  was  used  and  the  plant  was  changed  to  an  S-Function  implementing  the  model 
of  the  hexapod  according  to  Equations  (3.28)  and  (3.30),  on  pg.  36.  Simulations  were  made  to  excite 
the  piston  (bouncing),  tilt/tip,  twist  and  shear  modes  (see  Eigure  4.12). 

The  first  case  simulated  shows  the  performance  of  the  controller  for  the  vibration  axis  co¬ 
inciding  with  the  z-axis  and  no  rotational  component  (piston  mode),  with  a  frequency  of  40Hz-  The 
results  are  shown  in  Eigures  4.13  and  4.14.  As  one  can  see,  the  decentralized  controller  performed 
reasonably  well  and  no  competition  among  the  controllers  was  seen  in  steady-state.  The  nonlin¬ 
earities  in  the  plant  model  are  responsible  for  the  harmonics  seen  in  the  error  signal  observed  in 
Eigure  4.13(a). 
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Tilt  Twist 

Figure  4.12.  Hexapod  modes  [5]. 


The  seeond  ease  simulated  was  the  twist  mode  in  whieh  the  aetuators  needed  to  move 
symmetrieally  in  order  to  eaneel  the  vibration.  As  Figures  4.15  and  4.16  show,  the  performanee  was 
as  expeeted  and  was  the  same  as  that  obtained  in  the  piston  mode. 

In  the  next  ease,  a  more  general  vibration  disturbanee  was  simulated  by  making  eaeh  eom- 
ponent  of  the  signal  dn  have  a  different  amplitude  (namely  1,  3,  2,  —4,  3  and  —2m/s^).  Figures  4.17 
and  4.18  show  the  eontroller  performed  as  well  as  in  the  previous  two  cases. 
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The  last  simulation  added  a  seeond  tonal  {50Hz)  disturbanee  to  the  previous  ease  (ampli¬ 
tudes  2,1,6,— 12,1  and  —5m/s^).  This  disturbanee  was  not  assigned  to  the  eontroller  and  therefore, 
was  not  being  eontrolled.  The  amplitudes  were  the  same  as  those  for  Case  3.  It  ean  be  seen  in 
Figure  4.20  that  the  eontroller  performed  as  expeeted. 
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Figure  4.14.  PPH  Simulation,  Case  1,  Noisy  {SNRd  =  50dB)  (Accelerometer#!). 
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Figure  4.15.  PPH  Simulation,  Case  2,  No  Noise  (Accelerometer  #1). 
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Figure  4.16.  PPH  Simulation,  Case  2,  Noisy  {SNRd  =  50dB)  (Accelerometer  #1). 
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Figure  4.17.  PPH  Simulation,  Case  3,  No  Noise  (Accelerometer  #1). 


Frequency  {Hz) 


(a)  Error 


(b)  Weights 

Figure  4.18.  PPH  Simulation,  Case  3,  Noisy  {SNRd  =  50dB)  (Accelerometer  #1). 
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F.  COMMENTS 


The  simulations  verified  the  theoretieal  results  and  the  Adaptive  Disturbanee  Caneeller  ex¬ 
hibited  good  performanee  even  when  used  in  a  distributed  eonfiguration  and  in  the  presenee  of 
strong  nonlinearities.  It  is  also  important  to  mention  that  the  vibration  suppression  was  not  de¬ 
graded  by  the  presenee  of  noise  in  the  measurements.  This  result  was  eonsistent  for  both  the  SISO 
and  the  MIMO  simulations. 

The  simulations  also  eontirmed  the  theoretieal  results  about  the  unassigned  frequeneies:  the 
eontroller  does  not  affeet  the  error  at  frequeneies  different  from  the  assigned  ones,  indieating  that 
this  eontroller  is  very  seleetive. 

It  is  important  to  reeognize  that  the  properties  mentioned  in  the  last  two  paragraphs  are 
interrelated:  one  exists  beeause  of  the  other  and  viee-versa.  Sinee  the  eontroller  uses  a  synthetie 
signal  as  input,  it  is  very  insensitive  to  noise.  But  a  synthetie  referenee  signal  assigned  to  an  ineorreet 
frequeney  will  prevent  the  eontroller  from  eaneeling  the  disturbanee. 

As  expeeted,  the  eontroller  did  not  generate  any  new  frequeney  in  the  SISO  simulations, 
whieh  used  a  linear  plant.  The  Preeision  Pointing  Hexapod  model  used  in  the  MIMO  simulations 
eonsiders  the  geometrie  nonlinearities  of  the  hexapod.  These  nonlinearities  generated  harmonies 
that  appear  over  the  entire  range  of  frequeneies,  although  the  harmonies  are  very  small  beeause 
the  displaeement  of  the  top  plate  is  small  and  thus  the  angles  are  small.  Nonetheless,  this  effeet 
is  expeeted  to  be  amplified  if  other  nonlinearities  sourees  are  present,  sueh  as  statie  frietion  and 
baekslash. 

Note  that  one  important  limitation  of  this  method  is  the  requirement  for  the  preeise  fre¬ 
queney  of  the  disturbanee  or  the  performanee  will  suffer  eonsiderably.  This  effeet  was  already 
observed  by  Tsuei  et  al.  in  [33],  whieh  presented  a  modified  version  of  the  Filtered-X  LMS  with 
synthetie  input. 
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V.  COMPUTATIONAL  COST  ESTIMATION 


A.  OVERVIEW 

This  chapter  compares  the  computational  cost  required  to  implement  the  Adaptive  Distur¬ 
bance  Canceller  and  the  Multiple-Error  LMS  method.  Both  computational  complexity  and  memory 
requirements  are  compared. 

Section  B  defines  the  methodology  followed  to  estimate  the  computational  cost.  Section  C 
and  Section  D  present  the  computational  cost  for  the  Multiple-Error  EMS  and  the  Adaptive  Distur¬ 
bance  Canceller  respectively. 

Section  E  compares  the  computational  cost  and  memory  requirements  for  both  controllers. 
One  example  is  presented  in  which  a  hexapod  is  used  to  suppress  three  tones  with  one  harmonic 
each. 

B.  BASIC  CONCEPTS 

Vibration  control  requires  large  amounts  of  computational  power.  The  frequencies  involved 
are  usually  higher  than  the  traditional  control  methods;  therefore,  the  sampling  rates  are  also  higher. 
One  of  the  main  purposes  of  this  research  is  to  study  a  vibration  method  that  requires  minimal 
computational  resources.  The  controller,  after  being  implemented,  should  consume  the  least  amount 
of  electrical  power  possible.  Additionally,  the  computer  system  should  require  as  little  space  as 
possible,  since  electrical  power,  weight  and  space  are  very  expensive  in  space  applications.  Ideally, 
a  single  low-power  microcontroller,  with  no  external  memory  and  no  external  A/D  converters  should 
be  used. 
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Floating  point  operations  used  to  be  the  most  expensive  operations;  eonsequently,  the  eom- 


plexity  of  an  algorithm  is  traditionally  measured  by  eounting  floating-point  operations.  It  was  not 
uncommon  for  a  simple  floating-point  multiplication  to  be  emulated  in  software  or  to  need  a  ded¬ 
icated  external  chip  to  be  executed.  The  advances  in  microelectronics  changed  the  picture:  CPUs 
have  evolved.  The  new  DSP  CPUs  can  perform  one  floating-point  multiplication  as  fast  as  they  can 
perform  one  floating-point  summation  or  even  integer  math.  Table  5.1  shows  the  specifications  of 
an  embeddable  microcontroller/DSP  (Hitachi’s  SH7410,  a  CPU  for  embedded  devices).  This  mi¬ 
crocontroller  CPU  has  a  DSP  unit  embedded,  interrupt  controllers,  memory,  A/D  converters,  timers, 
etc.  The  power  consumption  is  very  small  (0.5W  versus  about  60W  for  an  Intel’s  Pentium  III  or 
an  AMD’s  Athlon).  Therefore,  if  a  vibration  control  solution  can  fit  in  a  few  dozen  kilobytes  of 
ROM  and  does  not  need  more  than  a  few  kilobytes  of  RAM,  the  savings  in  space,  weight  and  power 
are  very  significant.  As  an  even  simpler  alternative,  the  controller  could  be  converted  to  use  integer 
arithmetic  or  floating-point  emulation  and  use  a  CPU  that  does  not  support  floating-point.  These 
are  more  easily  available  and  come  with  a  wider  array  of  features  (some  featuring  eight  12bits  A/D 
conversors  and  eight  D/A  conversors),  making  it  possible  to  implement  the  controller  with  very  few 
components  and  little  power  consumption. 


Voltage 

3.3V 

Operating  current 

110  to  210mA 

Operating  power  consumption 

0.36  to  0.69W 

Standby  current 

SjuA 

Standby  power 

lOtiW 

Time  per  operation  (fetch,  multiply  and  divide) 

33ns  /  operation 

ROM 

A^Kb 

RAM 

SKb 

Cycle/instruction  basic  instructions 

1  cycle 

Cycles  needed  for  DSP  multiplication 

1  cycle 

Cycles  needed  for  integer  multiplication 
Interrupt  controller  embedded 

1  to  4  cycles 

Other  Features 

Timers,  A/D  converters,  serial 
communication 

Table  5.1.  Main  Characteristics  of  the  SH7410  (SH-DSP,  SH2  Family). 
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In  the  light  of  the  realities  of  the  new  CPUs,  evaluating  an  algorithm  purely  by  its  floating¬ 


point  eomplexity  is  elearly  not  enough  anymore.  Using  a  single  measure  is  no  longer  possible,  sinee 
the  performanee  of  an  algorithm  will  depend  heavily  on  the  implementation  and  on  the  arehiteeture 
of  the  hardware. 

The  goal  of  the  eomputational  eost  evaluation  in  this  researeh  is  more  ambitious  than  just 
eomparing  different  methods.  It  is  intended  to  be  used  as  a  guideline  to  evaluate  if  a  partieular  CPU 
ean  be  used  to  implement  the  algorithm  in  real  time.  As  sueh,  the  evaluation  of  eomplexity  estimate 
will  not  measure  only  floating-point  operations.  The  following  items  will  be  eounted: 


•  Floating-Point  Multiplication/Division  fMy  j  .  These  are  some  of  the  most  important  oper¬ 
ations  as  their  eost  ean  be  mueh  higher  than  most  other  operations  when  the  eomputer  does 
not  have  floating-point  operations  implemented  in  hardware. 

•  Floating-Point  Sum/Subtraction  fS f).  Important  for  the  same  reasons  as  presented  above. 

•  Integer  Multiplication  fM,- j.  Integer  operations  ean  be  very  fast  in  eomplex  CPUs.  But  the 
low-eost  and  low-power  ones  ean  usually  be  slower  for  multiplieation,  taking  several  eyeles. 

•  Integer  Division  fB,-).  Integer  divisions  are  a  eomplex  ease.  Some  low-end  CPUs  may  not 
implement  a  full  set  of  division  operations.  Thus  some  sort  of  software  emulation  may  be 
needed. 

•  Integer  Sum/Subtraction  (E>i).  These  operations  are  among  the  fastest  ones,  but  also  among 
the  most  frequent.  For  fast  CPUs,  they  are  as  fast  as  multiplieations  and  divisions.  But  for 
low-end  CPUs,  summations/subtraetions  are  mueh  faster. 

•  Matrix  Element  Access  (Matr)-  Almost  every  operation  in  a  eomputer  oeeurs  in  the  CPU,  and 
most  involve  registers  themselves  or  register  and  memory.  Fetehing  data  from  memory  to  reg¬ 
isters  takes  time,  whieh  is  very  important  for  matrix  manipulations.  Even  taking  into  aeeount 
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that  matrix  operations  frequently  use  the  same  base  address,  they  still  involve  mathematieal 
operations  on  the  index  prior  to  fetehing  the  aetual  value. 

•  Transcendental  Functions  (^ranc)-  The  aetual  implementation  of  these  funetions  varies  de¬ 
pending  on  the  algorithm  used.  Some  dedieated  hardware  are  ean  of  performing  some  of  these 
operations,  but  usually  they  require  several  CPU  eyeles.  While  this  may  not  be  pereeptible  in 
fast  parallelizing  CPUs,  the  eost  ean  be  quite  high  for  simpler  ones.  Operations  not  provided 
by  the  CPU  must  be  emulated  in  software  and,  again,  the  eomputation  time  is  very  dependent 
on  the  algorithm,  implementation  and  desired  preeision. 

C.  MULTIPLE-ERROR  LMS 

A  referenee  must  be  established  in  order  to  evaluate  the  numerieal  eomplexity  of  the  pro¬ 
posed  vibration  suppression  method.  The  most  natural  one  is  the  Multiple-Error  LMS  due  to  its 
widespread  aeeeptanee. 

The  ultimate  objeetive  of  any  vibration  eontrol  algorithm  is  to  generate  inputs  for  the  plant 

//,  whieh  in  the  Multiple-Error  LMS  is  done  by  using  the  Equation  (4.4),  reprodueed  below: 

e-i 

W  =  Yj  W  (5-  1) 

^=0 

where  x[n]  is  the  referenee  signal,  Wq^i  [n]  is  the  weight  matrix  and  Q  is  the  order  of  the  filter  and  n, 
is  the  input  of  the  plant  (0  <  /  <  A,  4- 1). 

The  number  of  operations  neeessary  is 

r  e-i 

=  N<  ^ 

I  ?=o 

=  (G- l)§/+2N{w^,,  [n]x[n-^]} 

=  {Q-\)Ef+Q{\Mf  +  mi+2Matr) 

=  (G-l)§/+GM/  +  G§;+2GML^r,  (5.2) 
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where  N{«}  means  number  of  operations  required  to  evaluate  •. 


Therefore,  the  number  of  operations  required  to  ealeulate  the  whole  veetor  u  [n]  is 


N{u[n\}=Ni{Q-l)Sf+NiQMf  +  NiQSi+2NiQMatr.  (5.3) 


The  next  step  is  to  ealeulate  the  weights  at  time  n+\.  To  do  that,  one  of  the  intermediate 


steps  is  to  evaluate  rg^i  for  all  q’s: 


N{ro,i[u]}  =  n|  ^ 

=  PS/+ (P+ l)N{/io,,[p]x[u-p]} 


=  PSf+{P+l){lMf  +  lSi+2Wl,,r) 


=  PSf+{P+\)Mf  +  {P+\)Si  +  2{P+\)Wi^tr,  (5.4) 


where  P  is  the  order  of  the  estimate  of  the  plant  and  ho^i  [n]  is  the  estimate  of  the  o*^  output  of  the 
plant,  due  to  the  input. 


Therefore,  the  weight  update  operations  eount  is 


r  K-\  ] 

N { Wq^i  [n+\]}  =  N <  Wqq  [n]  +2p  Y,  [n]  roj  [«  -  ^]  > 

=  IM^fr  +  1S/+  n|  2q  Y  [n]  roq  [«  -  ^]  | 

r  K-i  ] 

=  IMaf;.  +  1S/+  IM/  +  N <  Y  ^o[n]  ro,i  [n-q]  > 

=  IMflf;.  +  1S/+  IMy  +  {No  —  1)  §/+ A^oN{eo  [n]  rgq  \n  —  cf\} 


—  IM^f^  +  1§/+  IMy  +  {No  —  1)  §/  + 


(IMy  +  !§,•+  lMa,r  +  N{ro,,  [n  —  ^]}) .  (5.5) 


The  eomputational  effort  to  evaluate  roq[n]  is  given  by  Equation  (5.4)  and  is  P^f  + 
{P+  1) My  +  (E+  1) §,  +  2{P+  l)Ma(r.  Sinee  Equation  (5.5)  needs  time  to  be  shifted  to  u  — 
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one  integer  sum  must  be  added  to  the  result 


1]}  —  +  1S/+ IMy  +  (A^O  —  1)  §y  + 

A^„My  +  No^  i  +  N o'Mat  r  + 

No{P^f+{P+l)Mf  +  {P+\)^i+2{P+\)Matr  +  ^r) 

=  (iVo  +  1  +  2No  (P+  1))  ML,,  +  (1  +  {No  -  1)  +  NoP)  § /+ 

{\+No  +  No{P+l))Mf  +  {No  +  No{P+l)+No)^i 

=  (l+iVo(l  +  2(P+l)))ML„+A^o(/’+l)§/+ 

(l+iVo(/’  +  2))MI/+iV„(P  +  2)§,-.  (5.6) 

Therefore,  it  is  possible  to  evaluate  the  eomputational  effort  to  evaluate  the  matrix  W  (of 
dimensions  QxNi)  by  using  Equation  (5.6) 

N{W[n+l]}  =  iV,(2(l+^o(l  +  2(E+l)))ML„  +  A^,(2A^o(E+l)§/+ 

NQ  ( 1  +  (E  +  2) )  MI/  +  NiQNo  {P  +  2)  (5.7) 

Therefore,  the  total  eomputational  effort  required  to  implement  the  Multiple-Error  EMS 
algorithm  is  the  sum  of  the  effort  to  eompute  u  [n]  and  W  [n  +  1] 

N{W[n+l]}  +  N{u[n]}  =  iV;(2  (1  +  (1  +  2(E+ 1)))  ML„  +  (E+ 1)  §/+ 

( 1  +  (T’  +  2) )  MI/  +  NiQNo  (E  +  3)  § ,  + 

iV,  (e  -  1)  § /+ A^,'<2MI/  +  NQSi  +  2NiQMatr 
=  iV,(2(3  +  iV,(l  +  2(E+l)))ML„  + 
Ni{Q{l+No{P+\))-l)Sf+ 

NiQ{2  +  No{P  +  2))Mf  + 

NiQ{\+No{P  +  2))§,  (5.8) 
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Note  that  it  is  safe  to  assume  that  P  >>  1  and  2  >>  1  for  eomplex  plants.  Applying  these 


assumptions,  an  approximate  estimate  of  the  eomputational  effort  is 

f^{ME-LMS}  ^  A,-e(3  +  No(l  +  2(P+l)))M«f,+ 
Ni{Q{\+No{P+l))-l)Sf+ 

NiQ{2  +  No{P  +  2))Mf  + 

NiQ{l+No{P  +  2))§i 
=  2NiNoPQMatr  +  NiNoPQS  f + 

NiNoPQMf+NiNoPQSi 

=  NiNoPQ{2m„r  +  Sf+Mf  +  Ei).  (5.9) 

The  implementation  of  the  Multiple-Error  LMS  used  for  the  experiments  deseribed  in  See- 
tion  VI.D  uses  a  plant  model  of  order  30  and  a  filter  order  of  10.  Therefore,  an  estimate  of  the 
number  of  operations  would  be 

n{ME-LMS}  ra  NiNoPQ{2Wl,tr  +  ^f+Mf  +  E,i) 

ra  10,800(2M«f,  +  §/+M/  +  §,).  (5.10) 

Not  desiring  to  enter  into  the  details  of  the  CPU  used  in  the  experiment,  it  is  assumed  that 
the  CPU  takes  only  one  eyele  to  perform  eaeh  of  the  eomputed  measures.  Therefore,  the  total 
number  of  CPU  eyeles  needed  by  the  algorithm  for  eaeh  iteration  is 

n{ME-  LMS}  ra  NiNoQ^  (2  +  1  +  1  +  1)  =  54, 000.  (5.11) 

Sinee  the  sampling  rate  used  by  the  algorithms  in  the  Ultra-Quiet  Platform  are  is  IKEIz, 
that  means  that  the  CPU  would  need  to  exeeute  more  than  5410^  instruetions  per  seeond,  whieh 
gives  an  estimate  of  the  eloek  needed.  The  aetual  speed  needed  will  also  depend  on  the  eompiler 
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optimization,  on  the  code  implementation  and  on  other  activities  that  the  processor  may  be  executing 


in  parallel  (I/O,  D/A  conversions,  A/D  conversions,  filtering,  context  switching,  etc). 

The  microcontroller  described  in  the  previous  section  can  execute  each  instruction  in  33ns, 
or  about  3  10^  instructions  per  second.  Therefore,  this  microcontroller  cannot  be  used  to  implement 
the  Multiple-Error  LMS  on  a  hexapod  if  E  >  30  and  2  >  10. 

1.  Memory  Requirements 

An  estimate  of  the  algorithm  memory  requirements  is  very  important  to  decide  wether  a 
particular  hardware  setup  can  be  used  for  implementation  or  not.  For  microcontrollers,  this  estimate 
indicates  if  the  controller  will  fit  in  the  on-chip  memory  or  not.  For  more  complex  CPUs,  this 
quantity  is  important  to  verify  if  the  application  can  fit  in  the  cache  memory  (much  faster  than 
conventional  memory). 

From  previous  developments  it  is  clear  that  the  Multiple-Error  FMS  cannot  be  implemented 
in  a  low-power  CPU.  Since  a  full-fledged  processor  is  needed,  one  wifh  a  large  on-chip  cache  could 
be  selecfed  and  ideally  all  fhe  code  and  dafa  would  run  from  fhe  cache  and  fhe  speed  would  be 
improved  by  a  large  faclor. 

The  biggesf  portions  of  RAM  memory  required  by  fhe  Mulfiple-Error  FMS  mefhod  are 
used  fo  hold  fhe  plan!  Iransfer  funclion  esfimafe  H  and  fhe  weighl  W.  Thus,  fhe  number  of  values 
fo  be  sfored  is 


memh  =  {P+l)NiNo,  (5.12) 

as  fhe  plan!  has  order  P  and  dimensions  A,  x  No- 
Similarly,  fhe  weighl  malrix  W  demands 

memw  =  QNi-  (5.13) 
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The  total  memory  required  for  a  Multiple-Error  LMS  eontroller  is  approximately  the  sum 


of  the  memory  needed  for  the  estimate  of  the  plant  and  weights 

raerciMELMS  ~  Ni{Q+{P+  l)A/o)-  (5.14) 

The  number  of  elements  to  be  stored  for  the  hexapod  example  is  approximately 

mem?3  6x  (10+(30+  1)  x6)  =  1,476.  (5.15) 

The  total  memory  is  about  6^ytes  (4  bytes/float).  This  is  not  a  very  large  amount  and 
might  possibly  fit  in  the  eaehe  of  the  CPU  (if  one  is  present).  If  the  eode  size  ean  be  kept  to  a  small 
reasonable  amount  (a  few  kilobytes),  then  the  whole  eode  and  data  eould  run  off  the  eaehe  and  the 
speed  gains  would  be  substantial. 

D.  ADAPTIVE  DISTURBANCE  CANCELLER 

Cost  estimation  for  the  Adaptive  Disturbanee  Caneeller  is  mueh  easier  than  that  of  the 
Multiple-Error  EMS  approaeh.  As  in  the  previous  seetion,  the  eomputation  starts  by  evaluating  the 
input  to  the  plant  using 

N{M[n]}  =  N{a[n]eos(cOc?^)  sin(cOc?^)} 

=  N{a[n]eos(c0c«)}  +  N{6[n]sin(c0c?^)}  +  1§/— —  IMy 

=  2  (2My  -|-  T.ranc  +  )  -|-  ISy  —  ~  IMy 

=  3My  -|-  2T.ranc  +  +  l§/>  (5.16) 

where  a  [n\  and  b  [n]  are  the  weights  of  the  eontroller. 

Please  note  that  the  time  needed  to  feteh  the  variables  to  the  memory  was  added  as  a  matrix 
operation,  whieh  is  logieal  for  multiple-frequeney  and  distributed  implementations. 
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The  weight  update  equation  eost  ean  be  eomputed  as 


N{a[u+1]}  +  N{ft[u  +  1]} 


N{a  [u]  +qe[u]  eos((Orfu)}  + 

N{ ^  [u]  +  qe  [n\  sin(cOrfn) }  —  (eommon  part) 

2N { a  [u]  +  qe  [n]  cos((0^n) }  —  AWan  —  2My 
— 2Mj  —  IMatr  +  2§/  +  2N{qe[u]  eos(cOdu)} 

—  2Mj  —  IW^atr  +  2§y  +  2  (SMy  +  ITranc  + 

4My  +  2T  ranc  +  +  2E>f.  (5.17) 


By  summing  values  from  Equation  (5.16)  and  Equation  (5.17),  the  eost  is 


N{ADC}  —  My  +  2T^a„c  +  6MQf^  +  2§y+ 3My  +  2Tra„c  +  +  ISy 

=  4My  +  4Tranc  +  lOM^^r  +  3§/.  (5. 18) 

Transeendental  funetions  are  very  eostly,  and  it  is  reasonable  to  eompute  and  store  them  for 
use  in  both  parts  of  the  algorithm.  Ideally,  these  values  may  be  stored  in  registers.  However,  the 
worst  ease,  in  whieh  the  values  are  stored  in  eonventional  memory,  was  assumed.  Thus, 


N{ADC}  =  4M/  +  2T™„,+  14ML,,  +  3§/.  (5.19) 

Note  that  N {ADC}  may  be  halved  when  the  proeessor  used  to  implement  the  algorithm  has 
a  large  number  of  registers  or  temporary  variables  are  earefully  arranged. 

Now,  the  number  of  eomputations  required,  assuming  Nc  eontrollers  running  in  parallel  and 
Nf  frequeneies  to  be  eontrolled,  is 

N{ADC}^^^^^  =  iV,iV/(4My  +  2T™„,+  14ML,,  +  3§/).  (5.20) 

It  is  important  to  realize  that  the  above  estimate  did  not  assume  the  eode  to  be  optimized, 
and  the  aetual  eost  estimate  ean  be  improved  by  saving  temporary  results. 
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The  actual  count  obtained  for  a  hexapod  case  controlling  controlling  six  frequencies  (three 


fundamentals  and  three  harmonics)  is: 

N{ADC}6^6  =  36(4My+2T™„,+  14M^,,  +  3§/).  (5.21) 

Sine/cosine  functions  can  be  evaluated  by  using  a  series  expansion  or  a  table  look-up  if  the 
CPU  does  not  have  them  implemented  in  hardware.  The  series  expansion  would  take  around  300 
cycles  (?3  20  terms)  if  each  operation  could  be  performed  in  one  cycle.  Since  the  implementation 
only  needs  12  such  expansions  per  cycle,  the  number  of  operations  would  be 

N{ADC}6  6  36x  (4+14  +  3)  +  300x  12 

4,356.  (5.22) 

Thus,  the  total  time  would  be  approximately  when  the  CPU  mentioned  in  Table  5.1 
(p.  76)  is  used.  Note  that  the  calculation  of  the  transcendental  function  evaluation  is  responsible  for 
most  of  this  time  (82%).  A  table  look-up  would  easily  reduce  time  to  a  much  more  reasonable  value 
(less  than  lOOq^)  at  the  cost  of  added  memory. 

1.  Memory  Requirements 

The  Adaptive  Disturbance  Canceller  controller  requires  very  little  memory,  as  only  six  vari¬ 
ables  need  to  be  stored  and  some  of  these  can  be  shared  among  different  controllers.  The  total 
number  of  variables  needed  is 


memA?,.,]v^  =  4A,-A/  -f  A/  -f  1 .  (5.23) 

Therefore,  the  previous  example  would  need  151  floating-point  values,  or  about  604  bytes. 
Almost  any  microcontroller  that  is  fast  enough  to  implement  the  Adaptive  Disturbance  Canceller 
will  have  that  amount  of  RAM  embedded  in  the  chip. 
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E.  COMPARISON  AND  COMMENTS 


In  some  ways,  comparing  the  computational  and  memory  requirements  of  both  algorithms 
is  not  reasonable,  since  both  start  from  different  hypotheses.  The  Multiple-Error  LMS  approach  re¬ 
quires  both  a  reference  signal  and  a  complete  model  of  the  plant  transfer  function  whereas  the  Adap¬ 
tive  Disturbance  Canceller  requires  the  knowledge  of  the  disturbance  frequency,  no  reference  signal 
and  an  approximate  model  of  the  transfer  function  only  at  the  frequencies  of  interest.  Although 
there  is  some  overlap,  each  method  has  different  strengths.  In  addition,  note  that  the  Multiple-Error 
EMS  method  most  certainly  needs  a  full  computing  system  to  be  implemented  in  real  time,  which 
means  adding  memory  chips,  interrupt  controllers,  A/D  and  D/A  converters,  timers,  considerable 
space  and  much  electrical  power.  Eurthermore,  the  cost  increases  linearly  with  the  plant  model 
and  the  order  of  the  filter,  and,  as  a  result,  complex  systems  require  much  more  computations  than 
simpler  ones. 

On  the  other  extreme  is  the  Adaptive  Disturbance  Canceller,  for  which,  in  the  studied  exam¬ 
ple,  only  4,500  operations  per  iteration  (726  operations  per  assigned  frequency)  and  151  variables 
are  needed  for  the  six  actuators  controlling  six  frequencies  each  with  emulated  sinicos  calculations. 
The  Adaptive  Disturbance  Canceller  is  also  very  scalable:  the  complexity  increases  linearly  with 
the  number  of  frequencies  to  be  controlled,  but  remains  constant  with  the  plant  order. 

It  is  worth  mentioning  here  that  a  full  Adaptive  Disturbance  Canceller  solution  could  be 
implemented  for  a  hexapod  with  only  a  few  chips:  a  microcontroller  (which  would  include  timers, 
A/D  converters,  D/A  converters,  watchdog,  interrupt  manager,  serial  communication,  etc)  and  ana¬ 
log  anti-aliasing  filters.  Power  consumption  would  be  quite  low:  most  likely  below  IW.  Einally,  the 
whole  controller  would  be  small  enough  to  be  mounted  in  the  empty  space  between  the  actuators. 
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Quantity 

Expression 

Comp.  Cost 

ADC 

NcNf  (4My  +  ITranc  +  14Ma<r  +  3§/) 

ME-EMS 

NiNoPQ  {IMatr  +  §/+  My  +  §;) 

Mem.  Req. 

ADC 

4NiNf  +  Nf+\ 

ME-EMS 

Ni{Q+{P+mo) 

Table  5.2.  Computational  Cost  and  Memory  Requirement  Estimates. 


Quantity 

Adaptive  Disturbanee  Caneeller 

Multiple-Error  EMS 

Cost/iteration 

4,356 

54,000 

Memory  estimate 

604  bytes 

>  6,000  bytes 

Number  of  Frequencies:  three,  one  harmonic  each 


Number  of  Inputs:  six 
Number  of  Outputs:  six 
One  operation  /  cycle 
Floating-points  implemented  as  floats 
Table  5.3.  Adaptive  Disturbanee  Caneeller  x  Multiple-Error  EMS. 


87 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


88 


VI.  EXPERIMENTAL  RESULTS 


A.  OVERVIEW 

This  chapter  describes  the  implementation  of  the  Adaptive  Disturbance  Canceller  on  the 
two  hexapods  described  in  Chapter  II  (p.  9). 

Section  B  gives  an  overview  of  some  of  the  challenges  faced  when  solving  the  vibration  iso¬ 
lation  problem  using  the  Precision  Pointing  Hexapod.  These  challenges  ranged  from  inappropriate 
design  to  hardware  failure. 

Section  C  described  the  results  on  the  Precision  Pointing  Hexapod,  which  exhibits  strong 
nonlinearities.  Experiments  show  that  the  Adaptive  Disturbance  Canceller  algorithm  and  the  hexa¬ 
pod  configuration  are  not  sensitive  to  the  vibration  axis.  This  indicates  that  the  hexapod  and  the 
Adaptive  Disturbance  Canceller  are  adequate  for  a  generic,  off-the-shelf  vibration  isolation  mount. 

The  Ultra-Quiet  Platform  was  specifically  designed  for  vibrafion  isolafion  and  fhus  has 
heller  linearity  properlies  lhan  Ihe  Precision  Poinling  Hexapod,  making  if  a  heller  lesl  bed  for  com¬ 
parison  wilh  Ihe  Multiple-Error  EMS.  Section  D  (p.  104)  presenls  Ihe  resulls  of  Ihis  comparison. 

B.  EXPERIMENTAL  CHALLENGES 

This  research  was  lo  be  developed  entirety  on  Ihe  Precision  Pointing  Hexapod,  which  had 
nol  yel  been  delivered  when  Ihe  research  slarled.  Once  delivered,  several  lesls  were  conducted  lo 
verify  ils  behavior  and  also  lo  validate  Ihe  malhemalical  model  lhal  had  been  developed  prior  lo  Ihe 
availability  of  Ihe  hardware. 
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The  very  first  experiment  eonsisted  in  exeiting  the  bounee  mode,  in  whieh  all  the  aetua- 


tors  were  exeited  with  the  same  sine  wave.  The  output  of  eaeh  aeeelerometer  was  then  sent  to 
an  oseilloseope  and  to  a  speetrum  analyzer.  The  first  and  most  important  eonelusion  was  that  the 
system  exhibited  severe  nonlinearity.  Harmonies  up  to  2AKHz  were  observed  when  the  input  sig¬ 
nal  eontained  only  a  50Hz  sinusoidal  sampled  at  IKHz-  The  experiment  was  repeated  with  a  5Hz 
sinuwave  and  the  results  were  not  as  severe.  The  first  experiment  was  relevant  from  the  vibration- 
isolation  perspeetive  and  the  seeond  experiment  was  important  for  position  eontrol  (pointing).  The 
next  experiments  tested  the  shaker.  The  supplied  shaker  was  a  25W  eeramie  shaker  with  frequeney 
response  peak  at  42Hz-  The  useful  range  was  between  10  and  lOOHz- 

Realizing  potential  problems  with  aliasing,  the  first  experiments  were  made  using  lower 
frequeneies  (lOHz)  in  order  to  have  all  harmonies  below  500Hz-  Larger  travel  of  the  aetuators  was 
needed  at  these  lower  frequeneies  in  order  to  obtain  an  aeeeleration  large  enough  to  obtain  a  good 
SNR.  As  a  result,  the  shaker  required  more  power  and  operated  very  elosely  to  its  25W  speeifieation 
at  the  lower  and  of  the  frequeney  range  (10//z).  Exeiting  the  shaker  at  sueh  power  levels  inereased 
its  own  nonlinearity  effeets. 

The  design  of  the  hexapod,  as  delivered,  had  the  aeeelerometer  as  a  struetural  part  of  the 
strut,  as  shown  in  Figure  6.1.  In  this  design,  the  full  foree  aeting  on  the  strut  passed  through  the 
ease  of  the  aeeelerometer. 

While  running  the  first  experiments  (low  frequeney  due  to  aliasing  problems),  the  algorithm 
sometimes  beeame  unstable  and  the  aetuators  hit  their  limit  (40A  eaeh  aetuator).  Turning  off  the 
eontroller  manually  was  not  an  option  beeause  it  was  implemented  in  software  and  thus  it  failed 
when  the  software  entered  an  unresponsive  state.  Disabling  the  power  supply  of  the  aetuators  man¬ 
ually  was  not  easy  beeause  the  switeh  is  loeated  behind  the  power  amplifier.  Even  if  the  aeeess  were 
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Adapter 


easier,  it  still  depends  on  a  fast  response  by  the  operator.  Due  to  the  original  design,  the  foree  of 
the  impaet  was  transmitted  through  the  ease  of  the  aeeelerometers  and  two  of  them  ended  up  being 
damaged. 

Sinee  two  aeeelerometers  were  already  damaged,  all  experiments  were  halted  until  the  de¬ 
sign  flaw  was  eorreeted.  This  was  done  by  redesigning  the  adapters  so  that  the  aeeelerometers  were 
inline  with  the  strut  axis  but  no  foree  was  transmitted  though  its  ease.  The  new  adapter,  after  being 
assembled,  is  shown  in  Figure  2.7. 

Although  damaging  the  aeeelerometers  was  no  longer  a  eoneern,  preventing  damage  to 
other  eomponents  if  the  eontroller  beeame  unstable  was  still  a  eoneern.  The  range  of  the  aeeelerom¬ 
eters  was  ±2g  or  0.5  to  4.5V,  assuming  nominal  sensitivity.  It  was  then  assumed  that  if  the  voltage 
fell  outside  the  range  of  the  aeeelerometer,  then  an  unsafe  eondition  oeeurred  and  preventive  aetion 
had  to  be  taken.  The  aetion  seleeted  was  to  curtail  power  from  both  the  actuators  and  the  shaker. 
In  order  to  do  that,  a  circuit  was  designed  to  monitor  the  output  of  one  of  the  accelerometers.  A 
window  tested  the  voltage  and  cut  the  power  supply  of  the  actuators  and  the  shaker  if  the  voltage 
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was  beyond  safe  limits.  The  power  eould  only  be  restored  by  pressing  the  reset  switeh.  A  manual 


panic  button  was  also  ineluded  and  had  the  same  effeet  as  if  the  aeeelerometer’s  output  were  beyond 
aeeeptable  limits. 

Although  this  eireuitry  prevented  damage  to  the  experimental  hardware,  it  still  did  not 
solve  the  aliasing  problems.  Additional  hardware  filters  were  needed.  Among  the  several  options, 
switehed  eapaeitor  filters  were  seleeted.  These  filters  require  only  two  eapaeitors  and  one  inte¬ 
grated  eireuit.  The  cutoff  frequency  could  be  selected  by  an  external  oscillator.  Switched-capacitors 
anti-aliasing  filters  were  implemented  for  all  accelerometers. 

Once  the  safety  device  and  anti-aliasing  filters  were  implemented,  several  experiments  were 
conducted.  In  order  to  generate  acceleration  levels  high  enough  for  a  good  SNR,  the  shaker  had  to  be 
driven  very  closely  to  its  power  dissipation  limits.  To  solve  this  problem,  either  the  accelerometers 
or  the  shaker  had  to  be  changed.  Therefore,  the  shaker  was  replaced  by  a  larger  model  501T.  This 
reduced  both  the  overheating  and  the  shaker  nonlinearities. 

The  next  important  consideration  was  the  ball  joint  clearance.  The  manufacturer  stated  that 
the  clearance  was  around  2  —  4pm  for  each  joint.  From  the  simulations  and  experimental  results  it 
was  determined  that  the  amplitude  of  the  vibration  was  in  the  order  of  20  —  50pm.  Therefore,  the 
joint  clearance  was  found  to  be  the  next  problem  to  be  addressed.  The  best  solution  was  to  replace 
the  ball  joints  by  flextures.  Another  option  was  to  preload  the  joints. 

Since  preloading  the  joints  was  faster  and  less  expensive,  it  was  the  approach  first  tested. 
All  joints  were  preloaded  and  the  experiments  repeated.  Several  experiments  for  which  convergence 
was  not  achieved  previously  did  converge  after  preloading  the  joints.  This  was  the  final  setup  for 
the  experiments  conducted  in  this  research. 
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C.  PRECISION  POINTING  HEXAPOD 


The  hardware  setup  is  that  detailed  in  Seetion  II.C  (p.  15).  The  only  modifieation  was  the 
preloading  of  the  joints  with  springs,  in  order  to  reduce  the  backlash  present  in  the  joints  (around 
3/rm  per  joint,  according  to  CSA).  A  sinusoidal  signal  of  50Hz  was  sent  in  phase  and  with  the  same 
amplitude  to  all  actuators  to  evaluate  the  effect  of  this  backlash.  The  power  spectrum  density  of  the 
accelerometers  showed  lines  at  50Hz  intervals  until  around  lAKHz-  The  solution  was  to  preload 
the  joints  with  springs.  Unfortunately,  this  introduced  axial  forces  on  the  actuators  axis  which, 
according  to  Moltran,  created  more  static  friction. 

The  controller  box  was  implemented  exactly  as  defined  in  the  equations  given  in  Sec¬ 
tion  IV.C  (p.  46).  The  parameter  q  was  coarsely  tuned  for  each  frequency  being  controlled.  Adaptive 
Disturbance  Canceller  controllers  were  implemented  in  parallel  and  the  outputs  summed  in  order  to 
control  more  than  one  frequency. 

Several  different  modes  of  vibration  were  tested  in  order  to  evaluate  the  algorithm  sensitiv¬ 
ity.  The  definition  of  the  modes  is  shown  in  Figure  4.12  (pg.  63). 

1.  Along  Z-Axis 

In  this  experiment  the  shaker  was  mounted  directly  under  the  bottom  plate,  at  the  center. 
The  vibration  axis  was  coincident  with  the  z-axis. 

The  first  case  studied  was  the  single  tonal,  at  40Hz-  Harmonics  due  to  nonlinearities  were 
present  on  the  error  signal.  The  results  obtained  before  and  after  the  controller  are  shown  in  Fig¬ 
ure  6.2.  As  can  be  seen,  the  controller  did  suppress  the  controlled  frequency  {40Hz),  but  the  vi¬ 
bration  at  the  harmonics  (80  and  120Hz)  were  increased.  The  floor  noise  remained  unaltered,  as 
expected. 
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Figure  6.2.  Z-Axis,  One  Tone. 
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Figure  6.3  shows  the  results  when  40  and  80//z  frequeneies  were  eontrolled.  As  expeeted, 


the  eomponents  of  the  error  signal  at  both  frequeneies  were  redueed  to  the  noise  floor.  Note  that 
the  seeond  harmonie  {\20Hz)  inereased  up  to  about  30dB,  as  it  was  not  eontrolled.  This  effeet  was 
speeifieally  due  to  the  nonlinearities  of  the  hexapod. 

2.  Tilt-Tip 

Another  interesting  eonfiguration  was  the  tilt-tip  as  it  was  simple  to  understand  the  move¬ 
ment  and  the  displaeement  of  the  aetuators.  The  shaker  was  mounted  with  the  vibration  axis  parallel 
with  the  z-axis,  with  its  eenter  loeated  on  the  x-y  plane.  This  eonfiguration  did  not  involve  transla¬ 
tional  movement. 

Figure  6.4  shows  results  for  a  single  tone  in  whieh  the  seleeted  frequeney  {40Hz)  was  sup¬ 
pressed.  Note  that  the  statie  frietion  is  believed  to  be  responsible  for  the  inerease  on  the  harmonies. 
As  expeeted,  there  was  no  signifieant  ehange  in  the  noise  floor. 

Next,  the  system  was  setup  to  eontrol  disturbanees  at  39Hz  and  40Hz-  As  expeeted,  the 
eontroller  needed  a  longer  time  to  eonverge.  Another  important  finding  was  that  the  strong  statie 
frietion  exhibited  by  the  aetuators  generated  several  harmonies  multiple  of  IHz,  whieh  was  the 
differenee  between  the  two  frequeneies.  The  two  assigned  frequeneies  were  aetually  suppressed, 
but  the  energy  levels  on  the  uneontrolled  harmonies  were  inereased  (Figure  6.5). 

3.  Twist  and  Shear 

In  this  experiment  the  shaker  was  mounted  with  its  eenter  on  the  x-y  plane,  with  the  vi¬ 
bration  axis  perpendieular  to  it.  This  mainly  exeited  the  twist  mode.  The  results,  presented  in 
Figures  6.6  and  6.7,  were  as  expeeted. 
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Figure  6.3.  Z-Axis,  Two  Tones. 
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Figure  6.4.  TilFTip,  One  Tone. 
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Figure  6.5.  Tilt/Tip,  Two  Close  Tones. 
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Figure  6.6.  Twist,  One  Tone. 
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Figure  6.7.  Twist,  Two  Tones. 
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For  the  last  experiment,  shear  mode,  the  shaker  was  mounted  on  the  x-y  plane,  with  the 


vibration  axis  along  the  veetor  x  +  y,  where  x  and  y  are  defined  aeeording  to  Figure  2.1.  This 
induees  a  translational  movement  to  the  hexapod  (mueh  more  eomplex  than  a  translation  along  the 
z-axis).  The  results  shown  in  Figures  6.8  and  6.9  were  again  as  expeeted. 

4.  Comments 

It  is  very  important  that  the  algorithm  performs  reasonably  well  in  most  eases,  in  spite  of 
all  nonlinearities  exhibited  by  the  Preeision  Pointing  Hexapod.  It  is  also  important  to  mention  that 
the  sensors  used  in  this  hexapod  had  a  very  wide  range  ±2g,  and  thus  the  disturbanee  level  had  to 
be  quite  large  in  order  to  exeite  the  aeeelerometers. 

Another  very  interesting  eonelusion  is  that  the  hexapod  was  quite  insensitive  to  the  vibration 
type  (linear/angular)  and  direetion.  The  experiments  subjeeted  the  hexapod  to  disturbanees  from 
several  direetions.  The  performanee  of  the  eontroller  in  all  eases  showed  that  the  Stewart  Platform 
is  indeed  a  good  eonfiguration  for  a  generie  vibration-isolation  mount. 
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Figure  6.8.  Horizontal  Shear,  One  Tone. 
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Figure  6.9.  Horizontal  Shear,  Two  Tones. 
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D.  COMPARING  THE  PRECISION  POINTING  HEXAPOD  TO 


THE  MULTIPLE-ERROR  LMS 

It  is  very  important  to  compare  the  results  obtained  with  the  Adaptive  Disturbance  Canceller 
to  those  obtained  with  the  Multiple-Error  LMS  algorithm,  and  the  same  hardware  must  be  used  in 
order  to  perform  a  fair  comparison.  The  Ultra-Quiet  Platform,  described  in  Section  II.B  (p.  10),  was 
selected  because  it  exhibits  better  linearity  properties  than  the  Precision  Pointing  Hexapod. 

1.  Results 

First,  an  experiment  was  se  tup  with  a  single  tone  at  95Hz  in  order  to  verify  the  Ultra-Quiet 
Platform  characteristics  and  how  well  the  Adaptive  Disturbance  Canceller  performed.  The  results 
of  this  experiment  are  shown  in  Figure  6.10.  As  can  be  seen,  the  harmonics  were  amplified,  but  the 
desired  frequency  was  effectively  suppressed.  Figure  6.10(b)  shows  that  the  controller  did  not  affect 
other  frequencies  other  than  the  assigned  one. 

Next,  a  second  experiment  was  conducted  in  which  five  disfurbances  (65,  95,  125,  130  and 
\95Hz)  were  fo  be  confrolled  fo  verify  fhe  behavior  of  fhe  confroller  under  a  more  complex  sifuafion. 
LMS-based  fillers  do  nol  perform  well  under  Ihese  conditions  [15].  As  can  be  seen  in  Figure  6.11, 
fhe  Adaplive  Dislurbance  Canceller  was  able  fo  suppress  all  fhe  assigned  disfurbances  successfully. 
If  can  also  be  seen  lhal  fhe  error  signal  for  fhe  unconlrolled  case  shows  frequency  conlenls  olher  lhan 
fhe  five  tones.  These  are  believed  fo  be  a  consequence  of  fhe  cross-coupling  among  fhe  differenl 
frequencies  due  fo  nonlinearilies.  This  effecl  was  nol  observed  in  fhe  linear  simulations  and  was 
observed  in  fhe  simulalions  using  fhe  Precision  Poinling  Hexapod  model. 
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Figure  6.10.  Adaptive  Disturbance  Canceller  on  UQP  with  One  Tone. 
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Figure  6.11.  Adaptive  Disturbance  Canceller  on  UQP  with  Five  Tones. 

Having  verified  that  the  performance  of  the  Adaptive  Disturbance  Canceller  on  the  Ultra- 
Quiet  Platform  was  within  the  expected  parameters,  a  series  of  experiments  were  conducted  in  order 
to  compare  the  Adaptive  Disturbance  Canceller  method  with  the  Multiple-Error  LMS. 

One  important  characteristic  of  the  Multiple-Error  EMS  is  the  requirement  of  a  reference 
signal.  This  signal  is  used  both  to  generate  the  estimate  of  the  signal  r[n]  and  to  generate  the 
input  to  the  plant  through  the  filter  W.  It  is  expected,  therefore,  that  a  noisy  reference  will  degrade 
significantly  the  performance  of  EMS-based  methods.  With  this  information  in  mind  the  ME-EMS 
experiments  were  performed  with  both  noisy  and  noiseless  reference  signals  {SNR  between  21  and 
15dB). 

The  first  set  of  experiments  had  a  disturbance  with  a  single  tone  at  95Hz-  As  seen  in 

Eigure  6. 12,  both  the  Adaptive  Disturbance  Canceller  and  the  noiseless  Multiple-Error  EMS  had 
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almost  the  same  performanee  over  the  entire  speetrum.  It  was  expeeted  that  the  noisy  Multiple- 


Error  LMS  would  not  perform  as  well,  and  this  was  eonfirmed  in  the  experiment.  The  noise  floor 
was  signifieantly  raised  at  higher  frequeneies,  signaling  that  the  Multiple-Error  EMS  does  require  a 
elean  referenee  signal  in  order  to  perform  well. 

The  seeond  experiment  eonsidered  the  sum  of  two  elose  sinusoidals  (95  and  96Hz)-  Ade¬ 
quately  isolating  elosely-spaeed  tones  is  very  important  for  large  struetures.  The  beating  frequeney 
is  very  small  and  struetural  nonlinearities  will  demodulate  this  low-frequeney.  Sinee  most  struetural 
modes  are  at  low  frequeneies,  this  ean  eompromise  the  goal  of  the  vibration  isolation.  As  seen 
in  Eigure  6.13,  neither  of  the  methods  was  able  to  suppress  both  disturbanees  eompletely,  and  the 
Adaptive  Disturbanee  Caneeller  had  the  worst  performanee.  One  of  the  explanations  for  the  poorer 
performanee  of  the  Adaptive  Disturbanee  Caneeller  was  the  tuning  of  the  learning  eoeffieient  q.  The 
eoeffieient  was  eoarsely  tuned  and  was  not  optimized  for  any  partieular  disturbanee  eonfiguration. 

Another  important  faet  that  the  experiment  uneovered  was  that  the  floor  noise  of  the 
Multiple-Error  EMS  with  noisy  referenee  was  about  20  —  35dB  above  the  noise  floor  of  the  Adaptive 
Disturbanee  Caneeller  and  the  Multiple-Error  EMS  with  noiseless  referenee. 

The  last  experiment  set  eompared  both  the  Adaptive  Disturbanee  Caneeller  and  the  noisy 
Multiple-Error  EMS  when  eontrolling  three  frequeneies.  As  seen  before,  the  Adaptive  Disturbanee 
Caneeller  introdueed  several  peaks  in  the  power  speetrum  of  the  error  due  to  eross-eoupling  intro- 
dueed  by  nonlinearities.  Eigure  6.14  shows  that  the  Adaptive  Disturbanee  Caneeller  did  introduee 
several  terms  due  to  eross-eoupling,  some  about  l0  —  20dB  above  the  noise  floor.  The  noise  floor 
was  also  raised  by  5  —  lOdB  eompared  to  the  uneontrolled  error.  Although  both  methods  were  able 
to  suppress  the  disturbanees,  the  wide-hand  noise  introdueed  by  the  Multiple-Error  EMS  with  noisy 
referenee  ranged  from  20  to  40dB. 
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2. 


Comments 


As  seen  from  the  plots,  the  Adaptive  Disturbanee  Caneeller  performanee  is  very  similar 
to  the  performanee  of  the  Multiple-Error  LMS  with  noiseless  referenee.  When  noise  was  injeeted 
in  the  disturbanee-eorrelated  signal,  the  Multiple-Error  EMS  performanee  suffered  signifieantly: 
although  the  aetual  disturbanees  were  suppressed,  the  noise  floor  was  raised  eonsiderably.  The 
Adaptive  Disturbanee  Caneeller  was  not  affeeted  by  this  behavior  beeause  the  information  about  the 
disturbanee  (frequeney)  is  assumed  to  be  known  a  priori  and  thus  the  referenee  ean  be  synthetized 
without  noise  interferenee. 

Another  very  important  faet  was  cross-eoupling  due  to  nonlinearities  among  the  several 
frequeneies  verified  in  the  experiment  when  several  frequeneies  were  involved.  Although  this  ef- 
feet  was  first  notieed  using  the  Adaptive  Disturbanee  Caneeller,  it  also  oeeurred  in  the  noiseless 
Multiple-Error  EMS.  The  most  important  eonsequenee  is  the  inability  of  the  studied  methods  to 
deal  with  several  simultaneous  frequeneies  (seen  in  the  experiment  shown  in  Eigure  6.15).  This 
means  that  the  most  effeetive  vibration-isolation  approaeh  is  to  isolate  the  vibration  at  its  souree, 
where  fewer  frequeneies  are  present  and  thus  aetive  vibration  isolation  is  more  effeetive. 
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Figure  6.12.  Comparison  with  ME-LMS  —  Single  Tone. 
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Figure  6.13.  Comparison  with  ME-LMS  —  Two  Tones. 
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Figure  6.14.  Comparison  with  ME-LMS  —  Three  Tones. 
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Figure  6.15.  Nonlinearity  Effect  on  the  ME-LMS  —  Eive  Tones. 
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VII.  SUMMARY  AND  CONCLUSIONS 


The  new  generation  of  optieal  payloads  imposes  strieter  requirements  on  the  vibration  lev¬ 
els,  ehallenging  the  designers.  These  higher  standards  make  it  very  diffieult  to  meet  the  requirements 
using  only  passive  isolation  and  thus  aetive  solutions  have  been  pursued.  The  Stewart  Platform  is 
identified  in  the  literature  as  the  most  promising  system  for  loeal  vibration  isolation,  as  it  ean  gener¬ 
ate  forees  and  torques  in  any  direetion.  As  added  bonuses,  hexapods  exhibit  very  good  stiffness  and 
high  foree-to-weight  ratio.  Finally,  hexapods  ean  also  be  used  to  provide  fine  poinfing  eapabilify  fo 
fhe  payload  if  fhe  aefuafors  are  eapable  of  long  sfrokes. 

Several  feehniques  for  vibrafion  isolation  have  been  developed  sinee  1981  when  Widrow 
and  Burgess  presenfed  fhe  original  Filfered-X  LMS,  buf  fhe  researeh  of  vibrafion  isolation  using 
hexapods  inereased  appreeiably  only  after  fhe  firsl  half  of  fhe  1990’s  when  faster  DSP  ehips  beeame 
available.  Unforfunafely,  fhe  eondifions  in  spaee  are  nol  as  forgiving  as  on  fhe  ground:  weigh!,  spaee 
and  eleefrieal  power  are  sfriefly  limifed.  Therefore,  alfhough  several  vibrafion  isolation  feehniques 
are  available,  fhe  eompufafional  resourees  required  limif  fheir  employmenf  on  eurrenf  spaeeerafl 
designs. 


In  1998  Berfran  and  Monforo  presenfed  a  very  simple  algorifhm  eapable  of  performing  vi¬ 
bration  isolafion  on  SISO  sysfems.  Taking  advanfage  of  fhe  fael  fhaf  fhe  hexapods  used  in  vibrafion 
isolafion  use  smarf  sfrufs,  eaeh  sfruf  ean  be  viewed  as  a  SISO  system,  and  Iherefore  fhis  simple 
mefhod  ean  be  employed  on  hexapods.  The  algorifhm  is  inexpensive  enough  fo  be  implemented 
wifh  only  a  small  fraelion  of  fhe  eompufafional  requiremenfs  of  fhe  de  facto  sfandard  Multiple-Error 
LMS.  A  deeper  review  of  fhe  liferafure  revealed  a  laek  of  a  sfabilify  analysis  for  fhis  simpler  mefhod 
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that  takes  into  account  a  generic  plant.  An  analysis  of  the  influence  of  non-controlled  frequencies 


on  the  performance  of  the  algorithm  was  also  missing. 

It  was  not  known  initially  if  the  algorithm  would  work  properly  on  a  hexapod  due  to  the 
strong  coupling  among  the  struts.  Therefore,  it  was  important  to  simulate  the  hexapod  in  order 
to  evaluate  the  Adaptive  Disturbance  Canceller  in  a  controlled  environment.  Due  to  the  unique 
characteristics  of  the  Precision  Pointing  Hexapod,  an  appropriate  model  could  not  be  found  in  the 
research  and  therefore  a  suitable  model  was  derived.  The  model  takes  into  account  all  geometric 
nonlinearities  so  that  it  can  also  be  used  for  pointing. 

The  Multiple-Error  LMS  is  used  a  as  reference  throughout  the  work  due  to  its  acceptance 
and  good  performance.  The  Multiple-Error  EMS  algorithm  was  re-derived  and  the  resulting  equa¬ 
tions  used  for  the  computational  requirements  estimation.  So  that  the  Adaptive  Disturbance  Can¬ 
celler  could  be  applied  in  spacecraft,  a  rigorous  analysis  of  stability  was  provided  and  a  guideline 
for  the  learning  coefficient  q  selection  was  also  included. 

Results  show  the  method  to  be  very  selective  and  not  to  affect  the  error  at  frequencies 
unassigned  to  the  controller.  As  a  result,  the  actuators  do  not  need  to  be  oversized  to  prevent 
saturation,  as  is  the  case  with  the  Multiple-Error  EMS  and  other  methods  that  use  a  reference  signal. 

Simulations  on  a  SISO  plant  were  conducted  to  validate  mathematical  results  and  confirm 
the  algorithm  high  selectivity.  Simulations  were  also  conducted  on  the  Precision  Pointing  Hexapod 
model  to  verify  the  influence  of  the  coupling.  Results  show  that  the  coupling  did  not  compromise 
the  performance  of  the  controller  significantly. 

A  full  set  of  experiments  was  conducted  on  both  hexapods.  The  Precision  Pointing  Hexa¬ 
pod,  due  to  its  ball  joints,  exhibited  large  nonlinearities.  Despite  the  severe  nonlinearities,  the 
method  did  cancel  the  assigned  frequencies,  although  the  amplitude  of  the  unassigned  harmonics 
increased. 
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Experiments  on  the  Ultra-Quiet  Platform  were  required  to  eompare  the  results  against  the 


Multiple-Error  EMS  algorithm.  The  proposed  method  has  approximately  the  same  performanee 
of  the  Multiple-Error  EMS  algorithm  when  no  noise  is  present  on  the  referenee  signal.  However, 
Multiple-Error  EMS  performanee  was  signifieantly  degraded,  generating  wide-hand  noise  on  the 
error  signal  when  noise  is  present  on  the  referenee  signal.  The  Adaptive  Disturhanee  Caneeller  does 
not  suffer  from  this  problem  beeause  no  referenee  signal  is  used. 

The  Ultra-Quiet  Platform  was  built  for  vibration  isolation  and  elose  attention  was  paid  to 
linearity.  Even  so,  multiple-frequeney  experiments  showed  that  both  methods  suffered  from  the 
nonlinear  effeets,  generating  several  harmonies.  The  effeet  is  more  dramatie  in  the  Multiple-Error 
EMS  implementation  when  noise  is  present  in  the  referenee  signal,  inereasing  the  noise  floor  of  the 
error  signal  at  higher  frequeneies  Both  methods  were  able  to  eaneel  the  disturbanees  even  in  the 
presenee  of  nonlinearities. 

Einally,  eomputational  eosts  for  the  Multiple-Error  EMS  and  the  Adaptive  Disturhanee  Can¬ 
eeller  methods  were  eomputed  and  eompared.  Traditional  eost  analysis  uses  floating-point  opera¬ 
tions  eount,  but  reeent  advanees  in  mieroeleetronies  deereased  the  eost  of  these  operations  and  other 
operations  must  also  be  eounted.  Additionally,  memory  requirements  for  both  methods  were  evalu¬ 
ated,  providing  the  designer  with  some  valuable  information  when  seleeting  the  hardware  to  be  used 
for  implementation.  The  eomputational  eost  required  by  the  Adaptive  Disturhanee  Caneeller  ean  be 
smaller  than  one  tenth  of  the  Multiple-Error  EMS  requirements  and  similar  results  are  obtained  for 
the  memory  requirements. 
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A.  CONTRIBUTIONS 


One  of  the  main  eontributions  of  this  work  is  the  approaeh  from  the  eomputational  effi- 
eieney  perspeetive.  Although  mueh  researeh  has  been  done  on  vibration  isolation,  no  work  that 
foeused  on  eomputational  effieieney  was  available  previously. 

Researeh  showed  that  the  method  seleeted  laeked  a  stability  analysis  for  a  generie  plant, 
and  a  stability  analysis  was  derived.  In  addition,  a  guideline  for  the  adaptation  rate  /r  seleetion  was 
presented. 

In  addition,  this  work  provides  an  analysis  of  erosstalk  interferenee,  predieting  the  inliuenee 
of  uneontrolled  frequeneies,  whieh  has  not  been  addressed  by  the  researeh  eommunity  until  now. 

A  fully  nonlinear  state-spaee  model  suitable  for  vibration  and  position  simulation  of  large 
and  small  angles  was  developed  for  the  Preeision  Pointing  Hexapod.  This  model  shares  the  large 
angles  eharaeteristies  found  in  the  roboties  literature,  without  the  restrietion  of  small  aeeelerations 
and  is  suited  for  the  geometry  of  the  Preeision  Pointing  Hexapod.  The  models  of  the  aetuators 
aeeelerometers  were  ineluded  in  the  model,  allowing  it  to  be  used  as  a  tool  for  simulation  and 
eontrol. 

Even  though  simulations  were  used  in  order  to  better  understand  the  method,  all  the  theo- 
retieal  results  were  validated  experimentally.  The  experiments  on  the  Preeision  Pointing  Hexapod 
revealed  the  robustness  of  the  method  to  nonlinearities.  These  experiments  also  reinforeed  the  faet 
that  the  method  is  robust  to  sensor  noise,  eonsistently  aehieving  suppression  level  below  the  noise 
floor. 

A  performanee  eomparison  with  the  Multiple-Error  EMS  algorithm  was  eondueted,  inelud¬ 
ing  eases  with  noisy  referenee,  whieh  substantially  degraded  the  behavior.  It  was  experimentally 
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shown  that  the  Adaptive  Disturbanee  Caneeller  delivers  performanee  eomparable  to  the  Multiple- 


Error  LMS  in  the  optimal  ease  of  noiseless  referenee  signal.  If  the  Multiple-Error  EMS  uses  a  noisy 
referenee  signal,  its  performanee  is  degraded  and  the  Adaptive  Disturbanee  Caneeller  offers  better 
results. 

Einally,  the  eomputational  requirements  analysis  of  both  the  Multiple-Error  EMS  and  the 
Adaptive  Disturbanee  Caneeller  methods  were  ineluded.  This  analysis  revealed  that  the  simpler 
Adaptive  Disturbanee  Caneeller  ean  deliver  at  least  the  same  performanee  of  referenee  method 
using  1/10  or  less  of  the  eomputational  resourees  when  the  frequeney  information  is  available. 


B.  SUGGESTIONS  FOR  FUTURE  RESEARCH 

Improving  on  the  Adaptive  Disturbanee  Caneeller  is  not  as  easy  as  it  might  seem,  sinee 
performanee  improvement  must  be  earefully  balaneed  against  simplieity.  The  main  researeh  ob- 
jeetive  was  to  have  the  simplest  vibration  suppressor  possible,  and  improvements  might  easily  add 
eomplexity  that  will  defeat  the  main  goal. 

The  Adaptive  Disturbanee  Caneeller  has  two  main  weaknesses:  the  need  for  preeise  knowl¬ 
edge  of  the  frequeney  and  a  frequeney-dependent  learning  rate  (p).  The  first  issue  ean  or  eannot 
be  relevant  in  praetieal  situations.  Several  subsystems  provide  as  outputs  their  frequeney  or  operate 
at  a  fixed,  known,  frequeney;  ofhers  do  nof.  If  fhe  frequeney  is  nol  readily  available,  fhe  nexf  besf 
solufion  would  be  fo  measure  fhe  frequeney  direefly,  buf  fhis  is  nol  always  possible.  Again.  Iwo 
allernalives  ean  be  used  when  fhe  frequeney  informalion  is  nol  available  and  eannol  be  measured 
direefly:  1)  esfimafe  fhe  frequeney  or  2)  adapl  fhe  algorilhm  fo  Iraek  fhe  dislurbanee.  The  firsl 
oplion  requires  signiheanl  eompulafional  resourees  fo  be  implemenled  wilh  good  preeision  and  in- 
Iroduees  delay,  whieh  impairs  fhe  response  of  fhe  eonlroller  fo  ehanging  dislurbanees.  Eor  fhis  elass 
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of  problems,  using  an  algorithm  that  can  track  the  frequency  by  itself  is  highly  important.  This  is 


one  of  the  major  areas  for  improvement  in  this  method,  and  although  much  time  for  this  research 
was  dedicated  to  solving  this  particular  problem,  no  acceptable  solution  was  achieved. 

The  other  main  problem  is  the  selection  of  the  learning  rate.  The  algorithm,  as  presented  in 
this  research,  uses  a  different  learning  rate  for  each  frequency.  Although  the  actual  value  of  q  is  not 
critical,  it  is  more  interesting  to  use  a  single  value  of  q  for  all  frequencies.  In  order  to  do  that,  the 
transfer  function  of  the  plant  at  the  assigned  frequency  must  be  either  known  or  estimated  online. 
This  information  could  be  used  in  several  ways  to  improve  the  method’s  characteristics.  Assuming 
that  the  plant’s  characteristics  change  slowly,  this  estimation  does  not  need  to  use  much  computing 
power  requirements,  although  the  memory  requirements  will  certainly  increase  appreciably. 

Finally,  the  mathematical  model  could  also  be  improved.  The  use  of  quaternions  would 
make  its  computation  faster  and  improve  its  usability.  Static  friction  should  also  be  included  in  the 
model  of  the  accelerometers,  which  would  provide  more  realistic  simulations  at  lower  frequencies. 
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